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I .  INTRODUCTION 

Recent  developments  in  control  theory  and  practice 
coupled  with  digital  electronics  advances  have  produced 
increased  interest  in  digital  power  plant  modeling.  Although 
conventional  marine  boiler  modeling  is  not  a  new  concept,  it 
is  still  an  area  of  increased  interest  due  to  the  relative 
lack  of  analytical  data.  This  deficiency  is  due  primarily 
to  the  extreme  environment  that  must  be  endured  by  sampling 
equipment.  Advances  in  data  collecting  and  processing  equip¬ 
ment  have  made  boiler  data  collection  practically  feasible, 
and  in  fact ,  nuclear  steam  generator  development  depends 
heavily  on  actual  data.  However,  it  is  not  economically 
feasible  for  a  single  manufacturer  to  retrofit  modern  data 
collection  systems  on  a  standard  marine  boiler.  The  majority 
of  marine  boiler  models  available  are  not  attractive  from 
the  control  engineer's  viewpoint  because  of  one  or  more  of 
the  following: 

a.  The  model  is  over-simplified  to  the  point  of  being  a 
"teapot"  model. 

b.  The  model  is  completely  developed  and  written  in 
Laplace  transforms  and  the  state  space  equations  are 
too  difficult  to  extract. 

c.  The  model  is  over-complicated  to  the  point  of  being 
computationally  inefficient. 
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d.  The  causes  and  modeling  of  "shrink  and  swell" 
phenomena  are  ignored. 

This  model  attempts  to  compromise  between  simplification 
and  complication  while  including  a  theory  and  model  for  shrink 
and  swell.  The  model  follows  the  path  specified  by  Fini  [1] 
restated  below: 

a.  A  general  D-type  marine  boiler  model  is  prepared  in 
Continuous  System  Modeling  Program  (CSMP)  language. 

b.  A  general  FORTRAN  program  is  used  to  prepare  the 
initial  conditions  needed  by  the  CSMP  model. 

c.  The  initial  condition  program  depends  only  on  data 
easily  obtained  from  the  manufacturer's  technical 
manual  and  engineering  handbooks. 


II.  MODEL  CONSIDERATIONS 


A.  BOILER  TYPE 

The  D-type  boiler  was  chosen  for  modeling  because  of 
its  comparative  standardization  among  manufacturers  coupled 
with  its  popularity  in  naval  ship  power  plants.  This 
particular  boiler  is  one  installed  in  the  LHA-1  class  of 
United  States  naval  ships.  Since  D-type  boilers  operate 
under  a  wide  range  of  geometry  and  output  conditions  the 
initial  condition  and  model  programs  are  designed  to  operate 
over  similar  ranges. 

B .  MODEL  NOMENCLATURE 

The  model  nomenclature  differs  from  reference  [1]  because 
of  a  need  to  clarify  the  model  equations  and  afford  readability 
to  the  actual  computer  programs.  The  following  objectives 
were  considered  when  developing  the  nomenclature. 

a.  The  equations  must  be  easily  read,  necessitating 
minimum  variable-definition  referrals  by  the  reader. 

b.  The  notation  used  in  the  model  development  must  be  the 
same  as  that  used  in  the  computer  programs. 

The  resulting  notation  combines  a.  and  b.  and  consists 
of  frequently  used  notations  for  variables  followed  by  a 
subscript  to  indicate  state  point  location  and/or  time  using 
Figure  1. 

Table  1,  a  tabulation  of  general  rules  and  examples, 
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should  suffice  to  allow  the  reader  to  follow  the  model 
development . 

C.  EQUATION  DEVELOPMENT  FORMAT 

The  dynamic  model  and  initial  condition  programs  are 
developed  simultaneously  with  corresponding  equations 
written  in  appropriate  proximity  to  each  other.  A  notation 
of  [IC]  following  an  equation  indicates  that  it  is  an  initial 
condition  equation,  while  a  notation  of  [ME]  indicates  a 
model  equation.  In  some  instances  the  equation  applies  to 
both  programs,  an  indication  of  [IC,  ME]  is  used  to  indicate 
this . 

In  as  many  instances  as  possible  equations  are  developed 
directly  from  elementary  heat  transfer  and  fluid  dynamic 
principles.  This  stems  from  a  desire  to  encourage  more  re¬ 
search  in  this  particular  field  and  to  foster  a  better  under¬ 
standing  of  marine  boiler  dynamics  among  students  of  marine 
engineering.  The  model  was  developed  to  be  read  and  used  by 
students  with  a  minimum  of  additional  research. 

D.  ASSUMPTIONS 

Various  assumptions  are  made  to  facilitate  either  model 
simplification  or  program  tractability .  The  latter  are  not 
always  desirable  and  further  research  or  different  programming 
techniques  may  obviate  the  need  for  them.  The  following 
major  assumptions  are  made: 

1.  The  energy  transferred  to  the  generating  tubes  is 
distributed  uniformly,  i.e.  "even  heating"  exists. 


2.  There  is  a  thorough  mixing  of  the  fluid  in  the  steam 
drum. 

3.  Steam  generation  occurs  only  in  the  boiling  section 
of  tho  screen  and  main  bank  risers. 

4.  All  generating  tubes  and  downcomers  in  a  particular 
circuit  have  the  average  geometry. 

Uniform  energy  transfer  is  assumed  to  facilitate  model 
simplification  and  is  justifiable  since  the  generating 
banks  are  designed  to  optimize  uniform  heat  transfer.  Mixing 
of  the  fluid  in  the  drum  is  a  result  of  the  turbulence  in  the 
drum.  The  third  assumption  results  from  a  lack  of  reliable 
experimental  data  and  the  knowledge  that  D-type  boilers  are 
natural  circulation  toilers  relying  on  a  density  difference 
in  the  circulation  loop  to  promote  flow.  Steam  production 
outside  of  the  generating  banks  would  decrease  the  density 
difference,  thus  inhibiting  natural  circulation.  The  geometry 
assumption  is  necessitated  by  model  simplification. 
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III.  MODEL  DEVELOPMENT 


A .  GENERAL 

The  boiler  is  developed  from  control  volumes  for  which 
the  appropriate  equations  for  energy  transfer  and  storage 
rates  are  written.  The  initial  condition  equations  are 
generated  by  setting  the  rate  (d/dt)  terms  equal  to  zero  and 
solving  for  the  desired  variables. 

The  model  is  developed  with  three  different  elements: 

1.  Furnace-side  heat  transfer 

2.  Water-side  heat  transfer 

3.  Water-side  circulation 

B.  FURNACE-SIDE  HEAT  TRANSFER 

The  gas  flow  path  can  be  followed  on  Figure  1  where  it  is 
represented  by  the  dashed  line.  The  fuel  and  air  mixture 
enters  the  furnace  at  state  point  Q  where  radiation  heat 
transfer  occurs  to  the  screen  riser  tube  metal.  The  flue  gas 
then  leaves  the  furnace  and  exchanges  heat  via  convection  with 
the  superheater  (state  point  R-S),  main  bank  riser  tubes 
(state  points  S-T),  and  economizer  (state  points  T-U)  in  that 
order. 

At  state  point  Q  the  fuel-to-air  ratio  and  fuel  and  air 
temperatures  are  relatively  constant.  This  allows  the 
sensible  heat  of  the  fuel  and  air  to  be  lumped.  The  mass  flow 
rate  into  the  furnace  is  the  sum  of  the  air  flow  and  fuel  flow 


L3 


M3  <9  -  MFqq  t-  /vi A3  3  [»e] 

The  initial  values  of  fuel  and  air  flow  at  a  specified 
operating  point  are  known.  Therefore,  the  corresponding 
initial  condition  equation  is: 

MQQ0  =/vtR33^  +  MA  QQ(j6  ITC} 

The  mass  flow  rates  through  the  remaining  sections  of  the 
boiler  are  the  same. 


MQ3  -  -  rviss  -  rn-n-  =  MuLc 

-M SSJ*  =WTT^  =MuUj2f  ^ 

The  total  heat  supplied  to  the  boiler  can  be  written  as 
the  product  of  the  mass  flow  rate  of  fuel  and  the  "lumped" 
fuel  heating  value. 


QQ  -  MFQq  *  FttV/ 


At  a  specific  operating  point  the  heat  released  to  the  furnace 
is  available  from  the  technical  manual  as  is  the  furnace 
volume.  The  total  heat  input  at  steady  state  is  the  product 
of  these  two. 

QQ0  -  Gjwtf  *  FUSVOU.  \£(l] 


This  allows  the  computation  of  the  fuel  heating  value. 

FHV  •=  qo^/mfgjq^  \xe] 
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Furnace  to  Screen  Risers 


The  heat  transfer  to  the  furnace  screen  tubes  occurs 
primarily  through  radiation,  i.e. 

^  ~  -Tw) 

which  is  the  Stef  an -Boltzmann  law-*-where  a  is  the  Stefan- 
Boltzmann  constant  defined  as: 

6"  -  x  \c  T  Qru. 

Therefore  the  heat  transfer  to  the  furnace  screen  can  be 
written 

=  S16MA*  *  +  4.«S 

-  (w  +  4(^.0)  jtC  *  4.9})  \meJ 

where 

S  16m  A  A  =  O'  X.  k 

At  a  given  operating  point,  technical  manual  values  are 
specified  for  the  furnace  heat  absorption  rate  and  the  area 
of  the  radiant  heat  absorbing  surface  of  the  furnace  screen. 
The  product  of  these  two  yields  the  steady  state  heat  transfer 
rate  which  in  turn  can  be  used  to  compute  SIGMAA. 


This  equation  assumes  that  all  the  radiation  from  the  combus¬ 
tion  flame  strikes  the  screen  tubes,  and  both  the  flame  and 
the  tubes  behave  as  black  surfaces. 


L 
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=  GfFURAQ  #.  AfcAO 
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Siamm  =  QLfi  /((trr$  -h^.^)  *  *  4.0  - 
(tw0  •+  460. 0)  *  *  4.0) 

An  energy  balance  equation  for  the  mass  of  flue  gases  in  the 
furnace  region  can  be  written  in  the  form 

-oL  (rv\  Cp~r)  -  4b+al  energy  >s*fpl;ed  -  energy 
cH"  'frdnS'ferreci  4o  screen  4xJoes  — 

Gr\ero^j  gas  !e3^*ng  -furnace. 
rsQ.'on 

or  ^ 

d-  (w/vssqe.  *  OS(fL  *  =  qq  -ql- mrji 

*  cqe.  *  (n^-T/Vvie) 


Note:  TRR-TAMB  is  the  absolute  temperature.  At  steady  state 
the  rate  terms  are  equal  to  zero.  Values  of  total  heat 
supplied  to  the  boiler,  heat  transfer  rate  to  the  screen 
tubes,  total  flue  gas  flow  rate,  and  flue  gas  temperature 
are  available  from  previous  equations.  This  allows  the  specific 
heat  of  the  flue  gas  in  the  furnace  region  to  be  computed. 

-  (q<^  -  ql$£ )  /  (mrk^  ^  (tqr^  -iwe))  \xcT] 

During  a  transient,  the  density  of  the  flue  gas  in  the  furnace 
and  the  specific  heat  of  the  furnace  flue  gas  can  be  regarded 
as  constants  near  a  steady  state  operating  point.  The  dynamic 
equation  can  be  written  as 


OTRR.  =  (g»G) -ql-  MfcK.  *  Cqft.  *  (rRR.  -  TAMS'))  / 

(mASS^R.  ■JK.CC^')  tMe] 

where 

MASSQR.  =  feFLUE.  -51^  FurtVOU  \X(£\ 


4 


f 

i 


Thus , 


Ti^ie.  - 


OrRje. 


or  in  CSMP-III  language 


Tee.  -  iMrseu^TRR^  5  otrje.^ 


t«=3 


The  value  for  TRRO  is  available  from  the  technical  manual 
2 .  Flue  Gas  to  Superheater 


For  the  remaining  tube  banks  heat  transfer  occurs 
primarily  through  convection.  The  energy  given  up  by  the 
flue  gas  in  flowing  over  the  tube  surfaces  is  represented 
by  the  general  equation 

energy  =  mass  flow  *  specific 

given  up  rate  of  heat  of 

flue  gas  flue  gas 

*  change  in  temperature  of  the  gas 


e . 


1 


rv\  Cp 
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At  steady  state,  the  energy  transferred  to  the  superheater 
tubes  from  the  gas  must  equal  the  energy  transferred  from 
the  tubes  to  the  steam. 

Qz<f>  =  QH0  [xcT] 

The  steady  state  values  for  superheater  inlet  and  outlet 
steam  temperature  and  pressure  are  specified  by  the  technical 
manual  along  with  the  superheated  steam  mass  flow  rate.  This 
allows  the  direct  computation  of  Q40  and  Q30. 

Q40  -  %  (l4/UN)^  -  HMM0)  [ic3 

HMMO  is  the  enthalpy  of  saturated  steam  corresponding  to  drum 
pressure,  computed  with  curve  fitted  equations  [2],  HNNO  is 
available  directly  from  the  technical  manual.  The  specific 
heat  of  the  flue  gas  in  the  superheater  region  can  now  be 
computed . 

CftS  -  QZ0  /  7 fc  -TSS0V)  fccH 

Thus  the  dynamic  equation  for  heat  transferred  from  the  flue 
gas  to  the  tube  metal  is  written: 

Q3  ~  5ft  ox S  3ft  (trx-Tss) 

The  heat  transferred  to  the  tubes  via  convection  may  also  be 
represented  by  the  Grimson  correlation  [3]  which  states  that 
for  cross  flow  over  tubes; 
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Jad.  -  c.  f  'U.n  d  Pr  13 
k.$.  v  ^  ' 

where 

U.O©  -  Veloc.iy 
d  -  4“uJo€-'  c&i  awei'S/"'* 

(2,  -  Co*q.sta*vP 

Pr  =  Pra^cMi  nu-mloeo 

Since  h  can  also  be  written 

h  -  c| 

then  / 


also 

Uoo  —  _m_ 

leading  to  the  general  equation 
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Around  a  particular  steady  state  operating  point  the  physical 
properties  can  be  regarded  as  a  constant  and  the  physical 


dimensions  are  constant,  i.e. 


n  Pr'^  =  Const  a  nf 


which  implies 


For  flow  over  tube  banks  with  differing  rows  of  tubes  and 
differing  tube  geometries  the  constant  (n)  varies  little  and 
the  average  is  approximately  .6  [3].  Therefore,  in  the  case 
of  the  superheater  and  main  bank  riser  tubes,  the  correlation 


is  used  in  the  form 


q  =  C 


AT 


or  for  the  superheater 


Q3  -  KRS  *  MftR.  *  .(q  *,  (TR.S  -TWW) 

where  TRS  is  the  average  flue  gas  temperature 

TRS  *=>  (t**.  -^TSS*)  /Z-tf 

Since  the  steady  state  values  for  Q30,  MRRO,  TRRO ,  TSSO,  and 
TWO  are  given  by  the  technical  manual  the  constant  XRS  can 
be  calculated. 

KftS  -  Q  *  *  *(o  *  (trz&  -  feQ 
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To  provide  the  dynamic  solution  of  superheater  flue 
gas  outlet  temperature,  the  two  transient  equations  for  the 
flue  gas  to  superheater  tube  metal  heat  transfer  are  equated 
resulting  in 

Tfts-  TWW  —  Xz  crs  /krs. 

TdR.  -  TSS 

where 

TRS  -  (r&iz.  -V 

Solving  the  above  for  TSS  gives 

tss  -  (yrr.  ^twu;)/(phxi+1)  [me] 

whe  re 

P«+n  -  Z,  *  MrR/L  *  *  #H  #.  C*s /krs  I_me3 

3 .  Flue  Gas  to  Main  Bank  Risers 

The  equations  for  the  flue  gas  to  tube  metal  energy 
transfer  in  the  main  bank  are  derived  similarly.  For  the 
equation  involving  flue  gas  specific  heat; 

noss  #.  cst  (tss— ttt) 

The  Grimson  correlation  equation  for  the  main  bank  risers  is 

4^5" -  KST  &  MSS  *  fi.(o  ik  (TST  -TYv) 
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where 


Tsr  =  (rss  +ttt)  /z.0 

By  equating  these  two  equations  as  was  done  with  the  super¬ 
heater  the  flue  gas  temperature  leaving  the  main  bank  can 
be  written 

TTT  ^(TSS  *  TyV)  / -V-l) 

where 

phtz  -  z*  hss  *  *  0.4  ^  cst/kst  eJ 

The  steady  state  heat  transfer  is  computed  differently  than 
that  for  the  superheater.  Since  the  steady  state  heat  trans¬ 
fer  to  the  screen  bank  has  been  computed,  the  energy 
transferred  to  the  main  bank  is  the  difference  between  the 
total  energy  transferred  and  the  desuperheater  and  screen 
bank  energy  transfer.  At  steady  state  the  mass  flow  rate  of 
steam  out  of  the  boiler  is  equal  to  the  mass  flow  rate  of 
feed  in  the  boiler  and  the  total  energy  transferred  is  the 
mass  flow  rate  multiplied  by  the  enthalpy  change  between  the 
outgoing  steam  and  incoming  feedwater.  Thus,  the  steady 
state  energy  transfer  equation  for  the  main  bank  is 

-HeefS)-^0  -qwfi  Ire] 

Q990  is  computed  in  the  water  side  heat  transfer  section. 

Since  Q50  =  Q60  at  steady  state  the  heat  transfer 
coefficient  KST  and  the  specific  heat  CST  can  be  computed. 

KST  -  /  (WSS0  (-TS70  -TYYjjzO)  fai] 

CST  /  (fV\S$0  *  (rSS0  -TTT-0))  \Tc] 
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Flue  Gas  to  Economizer 


4 . 

The  economizer  tubes  differ  from  the  generating  and 
superheater  tubes  in  that  they  are  finned.  A  reasonable 
finned  tube  heat  transfer  correlation  was  developed  by 
Weirmann,  et  al .  [4] 

j  =  U  Pr  */3 

Cp  £»\a* 

G  =  maximum  mass  flow  rate 
max 

j  =  correlation  constant 


Since 


W 


specific  heat  of  flue  gas 


the  above 


correlation  can  be  written 

Q  -  j  A  A.T  Cg  €1  max, 

Pr^ 


G _  is  a  function  of  mass  flow  rate  into  the  tube  bank  and 

max 

tube  bank  geometry.  Similar  to  the  case  of  the  Grimson  coeffi¬ 
cient  used  for  the  superheater  and  main  bank  tubes  the  physical 
and  geometric  properties  are  considered  constant  leading  to 
the  use  of  the  correlation  in  the  form 


Therefore,  for  the  economizer  the  equation  is 


Q7  -  KXU.  *  n\tT  *  (tTu.  -tkx) 


where 

TTUL  -  TTT  +-TU.U. 

2.0 

Paralleling  the  superheater  and  main  bank,  the  heat  transfer 
may  also  be  written 

Q7  =  MTT  CTU.  =3^  (TTT-  [fv\£] 

The  steady  state  equation  development  follows  that  of  the 
superheater  and  main  bank  sections  and  is  not  repeated  here. 
The  equations  are: 

QS0  =  MftAjzS  *.  (HCe^-HAAjzO  .IbJ 
0(1$  -  taj 

CTU.  -  GafilCmrtp  *.(ttt <t>  - Tuu.jzO')  [re] 

Kru.  =  Q7^/(MTTjzS  5K.CTYu.j!5  -  TAX,fzS))  [rcQ 

TTU^5  =  (ttt0  yTu.u.0)  /z.  $  [re] 

Knowing  the  values  for  KTU  and  CTU  developed  by  the  initial 
condition  program  allows  the  solution  of  the  flue  gas  outlet 
temperature  in  the  transient  analysis.  Thus, 

ime]  tuu.  'Crrr  *.(phi3  -Lfl)+Z.0ik  tm)/ 

(PHX3  -t-1.0) 

ime]  PH-13  =  2.0^  CLTu/kTU. 
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C.  WATER  SIDE  HEAT  TRANSFER 
1 .  Risers 

The  heat  supplied  to  the  riser  banks  is  used  both  in 
the  boiling  and  nonboiling  length  of  the  tube;  however  the 
assumption  is  made  here  that  the  heat  transfer  is  all  under 
fully  developed  flow  boiling  conditions.  Since  the  nonboiling 
length  of  the  risers  is  relatively  small,  little  error  is 
produced.  For  fully  developed  nucleate  flow  boiling,  Levy 
[5]  suggests  the  equation 

A  Psat/3  AT^at  gm 

\.n*X  *  *0^  Sec- tt*' 

By  lumping  the  constant  area  and  denominator  the  riser  equations 
can  be  written  in  the  forms 

G\  -  C  *  PSAfT  ^  &  AXSfl-T 3 

Therefore,  the  two  riser  equations  are 

-  KV  *  PSAT #  (4.  (75/3. 0)  5K  (tW-  TSfVr)*  *.3.0  LmeJ 

and 

Q<c  -  K.Y  *  few  -TSVt)*  *3.0 

for  the  screen  and  main  bank  respectively.  A  simple  energy 
balance  on  the  tube  metal  yields  the  metal  temperatures,  i.e., 

^  ^  dT  - 

dt 
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The  equations  for  screen  and  main  bank  metal  temperatures  are: 


o TW  =  v  -K  cpv)  [ms] 

TY\/  -  iN>T£tf.L  j  OT/i/ ) 

OTVV  -  (qs*-  qgT)/ (MASSY  CP/)  j>E] 

tyy  -  ik>t^l.(tw0  borvy^ 


The  heat  transfer  coefficient  KY  can  be  computed 
directly  from  tube  geometrical  data  available  in  the 
technical  manual; 

KY  -  Aiee/we /V7S5LE-  <f>(o  Ere] 

Since  Q10  =  Q20  at  steady  state,  the  dynamic  equation  for 
Q20  can  be  solved  in  a  steady  state  situation  for  KV, 

tic]  kv-  qigS/(psA T**  61/3)#  z.&) 

Q60  was  computed  in  Section  III-B.3.  This  allows 
the  calculation  of  the  initial  main  bank  riser  tube  metal 
temperature . 

[  ic]  Tyy0  =  (Q(o0/(  KY  *,  Ps*r  **  fa  fa)))  *  %  fa  fa)  » 

In  the  case  of  the  screen  risers,  the  initial  tube 
metal  temperature  is  available  from  the  technical  manual  and 
the  initial  heat  transfer  rate  has  been  previously  calculated. 
This  allows  the  calculation  of  the  screen  riser  steamside 
heat  transfer  coefficient, 

I  IC]  kv  =  ql0/(RSAT  *  *  (rVt/0-T SA-r)**  3.0) 
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2 .  Desuperheater 

The  Dittus-Boelter  correlation  [6]  is  used  to  compute 
the  heat  transfer  rate  from  the  steam  to  the  desuperheater 
tubes.  This  correlation  in  its  basic  format  is 


€ 


K 


where 

Nu  =>  Nusselt  number 
Re  =>  Reynolds  number 
Pr  =>  Prandtl  number 
n  =>  .3  (for  cooling  fluid) 

When  the  appropriate  variables  and  constants  are  substituted 
for  the  dimensionless  numbers,  the  Dittus-Boelter  correlation 
can  be  written 

0,^  K.rV*A,T 


where  the  constant  K  is  defined  as: 

^  ^  ~H>ermaA  IqeaA*  ■VrarSHer  area 

•Vubcu  <Ai 


/ _ d _  yg 

\  TT^;  -Vute  dUa^eta/'-  \Ascosi4y  %  nr  h*bes  ) 


*  Pr*' 


and  AT  is  the 
Thus , 


log-mean-temperature  difference  (LMTD). 
for  the  desuperheater,  the  equation  is 
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-  KfJP  ^  ^  LMTOKJP 


KNP  can  be  computed  directly.  Thus, 


knj 9  -  ((.  faz  *.  THdoiOfj)  /dmp)  DIOP 

*  Viscoij  *  tJTiAfiOS))#.*  <=KAU3K*0.3 

*  a  Reads 


The  energy  given  up  by  the  steam  to  the  tubes  may 
also  be  written  in  terms  of  steam  specific  heat,  flow  rate, 
and  temperature  difference. 

~  (LNJP  -rpp) 

Similarly  since  the  heat  transfer  to  the  water  in 
the  drum  is  totally  by  convection; 

QSS  -  k*  #.(tss  -thh) 

For  initial  condition  calculations  THE  is  considered 
to  be  equal  to  TSAT.  Therefore 

ks  -  QW0 /Crzutf-T sat)  [x<J 


The  desuperheater  outlet  temperature  is  solved  by 
equating  the  two  dynamic  equations  for  Q9  and  solving  for 


the  desired  outlet  temperature. 

TPP  a  (Tb)h)  -  TZ-2’')/ (exP  (mP  /  (C.KiP  *  MWU**  0.Z  ))) 


+T %Tk 


ImeJ 
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The  instantaneous  desuperheater  tube  metal  tempera¬ 
ture  is  obtained  by  an  energy  balance  on  the  tube  metal 
similar  to  that  for  the  riser  tube  metal. 

’  %  - 

Solving  the  appropriate  equation  for  dT/dt  results  in 

0T££  =  OPs)  [we] 

-  liOTSfU.  DTft?  £v\ej 

The  specific  heat  of  the  steam  in  the  desuperheater 
(CNP)  may  be  calculated  using  the  known  steady  state  values 
for  desuperheated  steam  flow  rate  and  the  desuperheater  inlet 
and  outlet  temperatures  and  pressures. 

~  *.(hnjnj0- tipptf)  \jc] 

(Lkjp  -Tppp))  [rcQ 

With  the  value  of  Q90  computed  the  initial  log-mean- 
temperature  difference  can  be  calculated. 

VJWTDrOP  a/P  *  MK)h)0  *  # 

The  log-mean-temperature  difference  is  a  function  of 
iteam  temperature  in,  tube  metal  temperature,  and  steam 
temperature  out;  therefore  the  initial  tube  metal  temperature 
may  be  calculated. 


rv> ass  jfc  dp  K  dT~ 

dt 
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T <j>  =  (tkJk3^  -  EXflODS  fcTTRR^)/ (i.0  -EXpoDS}  fel] 

expoos-  exp((tkjnj^  -tpp$  VuvxToiJP)  trcTJ 


3 .  Economizer 

Paralleling  the  desuperheater  development  the  Dittus- 
Boelter  correlation  is  used  to  relate  the  heat  transfer  from 
the  tube  metal  to  the  feedwater.  Thus, 

tOa  -*•  ,fiZ3  Ri%n 

the  only  difference  being  that  n  is  now  .4  vice  .3  because  the 
fluid  is  now  being  heated.  The  appropriate  constants  are 
again  lumped  yielding 

<3?  -  KX  &  WAA  *.  X  0.?  ^  UNAT0A6 

An  energy  balance  on  the  water  gives  the  heat  absorbed  in  two 
other  forms  ; 

-  ro AA  *  (h«S6  -  HAa) 

Q?  ~  wAA  jK  CAS  *  ("nSG  -TAa)  \>e] 


At  the  specified  operating  point  the  inlet  and  outlet  condi¬ 
tions  of  the  economizer  as  well  as  the  mass  flow  rate  are 
given,  allowing  the  computation  of  the  steady  state  heat 
transfer . 

*  maa^zS  *  -haa^)  Ire] 
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The  heat  transfer  coefficient  KX  is  computed  directly. 

kx  -  (( . tfzz -5K.thcoiua)/da8)  m  (4>0/(pr*oA8*< 

VlSCOA  *  NJTuSEC))**  T  *  PRAA**^-4 

*  EC,  £rc] 


This  permits  the  solution  of  the  steady  state  log-mean-tem¬ 
perature  difference. 

UmTOfte  =  QS-JzS/tKX#  r^AAjzS  <Z>.?) 

Following  the  identical  path  delineated  for  the  desuperheater, 
the  tube  metal  initial  temperature  can  be  computed. 

TXx^  .=  (TAA0-T6 *  £XPoZc)/(\.$  -EXP^EC,)  [xcj 
BCPOEC,  -  EXP((T8S^-TAA0  l^c] 


An  energy  balance  on  the  economizer  tube  metal  yields 


its  instantaneous  temperature. 

DTXX  *=  (q?-C?7)/(^ASS)C  *  c,pa)  [>e] 
T*X  »  DTxx) 


By  equating  the  formulations  for  Q8  involving  LMTD 
and  specific  heat,  the  economizer  outlet  temperature  may  be 


calculated . 

T6e=  (t-AA  -TXX.)/(ex,P  (k.x/(C.A&*MAA*.*  $.2.))) 

+  tx.*  tA'eH 
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The  specific  heat  of  the  feedwater  in  the  economizer 


is  calculated  at  steady  state  using  the  computed  value  of 
steady  state  heat  transfer  and  the  given  mass  flow  rates  and 
inlet  and  outlet  temperatures. 

4 .  Superheater 

The  development  of  the  superheater  equations  follows 
that  of  the  desuperheater  and  economizer  and  for  that  reason 
the  development  will  not  be  repeated.  The  equations  are 
listed  below. 

-  0\*j  %  mnw  *  (rwrJ  -tmm') 

(hWk/0  -  \zt\ 

dntJ  -Q^\<p  /  c] 

The  initial  superheater  log-mean-temperature  difference 
is  computed  directly  allowing  the  subsequent  calculation  of 
the  tube-metal-to-steam  heat  transfer  coefficient  used  in  the 
Dittus-Boelter  equation. 

L»*\TDrv\KJ  -  (tkJnJ0-tw0  V (*1-06  ((.TvJW<fi-TfAiA0)  / 

(tv/w0-tnJn;0  )))  \xc\ 

KW  =  *  *,  0.?  *  uvtpcmkj)  Ifc] 


An  energy  balance  on  the  tube  metal  yields  the  instantaneous 
tube  metal  temperature,  the  initial  tube  metal  temperature 
being  available  from  the  technical  manual. 

DrWW  -  (q3-q4)  / ^  cpw) 

TVv/V\/  “  OTIa/U/) 

The  superheater  outlet  temperature  can  now  be  calculated. 

tnjnJ  *  -twia/)/(exp  (kw/(cmiJ% 

i-TWu/ 

D.  WATER-SIDE  CIRCULATION 
1 .  General 

The  water/steam  side  circulation  equations  are  by  far 
the  most  difficult  to  visualize  and  codify.  These  equations 
must  be  accurate  in  order  to  predict  phenomena  such  as  shrink 
and  swell  while  appropriate  assumptions  and  simplifications 
must  be  made  in  order  to  make  the  equations  tractable. 

As  can  be  seen  in  Figure  1 ,  the  feedwater  enters  the 
economizer  at  state  point  A,  passes  through  the  economizer 
to  the  steam  drum  where  it  becomes  part  of  the  water  volume. 
The  liquid  leaves  the  drum  via  the  downcomers,  both  main  bank 
and  screen  at  state  points  G  and  C  respectively.  The  main 
bank  downcomer  delivers  its  liquid  to  the  main  bank  risers 
via  the  water  drum,  and  the  main  bank  risers  deliver  the  fluid 
to  the  steam  separators.  The  flow  through  the  screen  risers 
is  the  same  with  the  exception  that  there  is  no  water  drum  in 
this  circuit.  The  steam  separators  separate  the  majority  of 
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the  water  from  the  steam  leaving  the  drum  at  state  point  M 
and  the  majority  of  the  steam  from  the  water  being  discharged 
back  into  the  drum  liquid.  The  steam,  with  a  negligible 
amount  of  water  passes  through  the  superheater  via  state 
point  M-N.  At  the  outlet  of  the  superheater  the  steam  leaving 
at  path  III  is  used  for  "main"  steam;  that  leaving  via  path 
II  travels  to  the  desuperheater  where  it  is  cooled,  then 
leaves  via  state  point  P  and  is  used  for  auxiliary  steam. 

The  water  leaving  the  separators  with  a  small  amount 
of  entrained  vapor  is  mixed  with  the  incoming  feed  water  and 
forms  a  "fc  iming"  vapor/liquid  mass  in  the  bottom  half  of 
the  drum.  The  liquid  from  this  mass  leaves  via  the  downcomers 
and  is  circulated  through  the  loop. 

2 .  Downcomer  Pressure  Drop 

A  closer  look  at  the  downcomer  flow  is  necessary  to 
justify  assumptions  that  are  made  in  the  development  of  the 
pressure-drop  equations. 

Circulation  ratio  is  defined  as  the  ratio  of  the  total 
weight  of  steam  liberated  to  the  drum  [7],  For  a  600  pound 
marine  boiler  this  circulation  ratio  is  on  the  order  of  20:1 
[7  and  8] .  This  implies  that  for  every  21  pounds  of  liquid 
flowing  down  the  downcomers,  20  pounds  of  it  has  already 
traveled  up  the  riser  and  is  at  saturation  temperature.  There¬ 
fore,  assuming  the  downcomers  are  perfectly  insulated,  it  is 
reasonable  to  assume  that  the  downcomer  liquid  is  at  or  very 
near  saturation  temperature. 


The  momentum  equation  for  the  downcomers  may  be 


written : 

pressure  at  top  -  pressure  at  bottom 
=  frictional  loss  -  gravitational  head 
+  entrance  loss  +  bend  loss 
+  exit  loss  +  inertia  force 

The  inertia  force  term  may  be  considered  negligible  [9] . 

This  is  the  quasistatic  approximation  which  basically  states 
that  pressure  waves  move  much  more  rapidly  through  the  system 
than  the  important  time  constants.  This  is  mathematically 
equivalent  to  the  elimination  of  a  large  negative  eigenvalue. 
This  quasistatic  approximation  is  only  good  for  on-line 
transients  and  does  not  apply  for  extremely  large  discon¬ 
tinuities  which  would  result  in  the  boiler  being  taken  out  of 
operation,  (e.g.,  a  ruptured  tube).  Therefore  the  downcomer 
momentum  equation  can  be  written: 

pressure  at  top  -  pressure  at  bottom 
=  (friction  factor  *  downcomer  length 
-r  downcomer  diameter  +  entrance  factor 

+  bend  factor  +  exit  factor) 

2 

*  (downcomer  mass  flow  rate)  f  (2  *  (downcomer 

2 

cross  sectional  area)  *  density  of  fluid 
in  downcomer  *  g  )  -  density  of  fluid  in 
downcomer  *  gravitational  acceleration  * 
height  of  downcomer  *  g 


In  the  model  notation 

PdPO  -  (fCO*.  L<Lo/ocD  -I-  EtOTlRCO  +  6&dOC£ 

-hEXireo)  *  M<LC  *.  *  2.0/ (2*0  *.*<!£**  2.0 
*.  R14CCD  *  6C.)-  lertOCD^  <5  *  SCO  /sc 

P6PK  *  LGtt/OfiH  -HEn^GB  +  8E>Jd£B 

+  EXirsH');^  MSS  *  *  2.0/(2.05«.  /\SH**Z.0 
*  RHOSB  ^  Sc)  “*  KBOGB  *  £  *  ^GB /SC 


for  the  screen  and  main  bank  downcomers  respectively.  FCD  and 
FGH  are  friction  factors  of  the  form 

4  =•  1  /(l.74  -2,  lo<5  O/ks)) 


where  R  is  the  pipe  radius  and  KS  is  the  relative  sand  rough¬ 
ness  [11]  . 

Fto  »  -2.fi  *  ALosigS(ks.<io))  [iCj«e] 

F6H  =  4.0/ (l/W  -Z-0  *  10  (kS<SH»  [k.,Me] 

3 .  Riser  Pressure  Drop 

The  momentum  equation  for  the  riser  boiling  section 
must  take  two  phase  flow  effects  into  account  because  the 
flow  in  the  boiling  section  is  not  homogeneous.  Thus,  there  is 
a  relative  velocity  between  the  liquid  and  vapor  phases  here. 
The  steam  separators  are  included  in  the  riser  length.  However, 
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the  effective  length  of  the  separators  and  thus  the  pressure 
drop  are  considered  negligible  because  of  a  general  lack  of 
information  concerning  these  items. 

The  pressure  drop  due  to  single  phase  or  homogeneous 
flow  in  the  nonboiling  riser  section  can  be  written 


AP  =  AP  +  AP  +  AP 

SPF  acceleration  friction  gravity 


^acceleraf.dn  ~~  ^  ^  ~~ J3.*) 

AfV*/,W  =.  a  y.  Ke..^kf *.  Gc,Jr  +  ) 

Experiments  conducted  by  Babcock  and  Wilcox  [10] 
indicate  that  these  homogeneous  flow  pressure  drops  may  be 
modified  to  give  the  appropriate  two  phase  flow  pressure  drops 
using  correction  factors  that  are  functions  of  outlet  quality 
and  operating  pressure.  The  two  phase  flow  pressure  drop  can 


be  written; 


AP 


TPF 


=  AP  , 

accelerat ion 

*  r 

accelerat ion 

+  AP 

*  r  _  . 

f rict ion 

friction 

+  AP 

*  r 

gravity 

gravity 

where  the  r  terms  are  two  phase  flow  correction  factors 
available  by  using  curve  fits  of  the  data  from  reference  [10], 
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The  form  of  these  r  terms  become: 

RGOWE  -  2.4/W4  *  XFF0  *  *  2.0  -  (o.^<o(p  *  XFF0  f  .477<<? 

26eA Vk  *  XU.^  *  *.  2.0  -  <o.  S0t>(c  X  Xll^  4.477&  txcj 

RAclE  -  ir.HSV^XFF^*  *  2.0  -H?.4W  ¥XFF0  &ej 

^cut  =  k:4s&4#xil0*  *  z.f  tis.4W  *xll0  -.00^07  yxc] 

fcFSlCE  - -34. 0?2X  *  xFh0  *  *  2.0  +  23- 7ifc>4  *  XFF0 

4- .  *734  ^ 

<?FKicK  -  -34. 0?^  *  Xu0  *  *  2.0  4*23.  7 M  *  XU0  0^1 

4-.S734 


Therefore,  in  the  boiling  region  of  the  riser 

^^■nPF  3  ^  ^Vd  2,  j*  {bo.'I.Aj  )cr%gR\^rv>  P{^-c 

4-  o£  igeil:^  reg.orQg.o  ryd\/ 

% 

The  total  pressure  drop  across  the  length  of  the  riser  is; 


AP  APSpp  +  APTPF 


In  model  notation  the  equations  are 
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I 


POPF=*  (?v\00*jfc  2.0  *  (jRUOee-I^WOOD))/ (<36  4r  l^HceE  * 
fcUcOO  *  AOE  *  *  2.0)  4-(M.0  *  foe.  *  lDE 
*  MOO  *  *  2.  0  )  /  (<66  X  DOE.  *  (RHOOD  +• 
RHoeE)  X  AD£  *  X  2.  0)  4  fe  X£OE  *  (rhcdd 

4-fcrtOEE))/  (66  X2.0)  4-  (moo*#  2.0  * 
£AC.Le)/(66  *  ^HoEE*  AEF  *  X  2.0)  4- 
(2.0  *  FEF  *  LEF*  flV>D  *  *  2.0  -&KFW?] ce)/ 

(66  *  OEF  *  KBCEE  *  AEF  *  X  2.0)  4- 
(6  X  S£F  *  £HoEE  X  €6£A\/<e)  /  66 


^(Wr  #  *2.0  *  (^w<-^hott))/6C^^hokx.)K 
RBoTT*  ATK*  *  **0)f(H.0£ 

*-4>  +  LTK)/ (66  *  D3X  *  ((2HClT3’4^HoKK)^ 
ATX  **2.0  *6  *^TK.*(^Hc3tr-t-^oK^)/ 

(66  *  2.0)  +-{Wr*  *  2.0  ^  KKClK)/ 

(66  *  KBOKK  *  AKU  *  *  2.0)  *  (2.0  * 
FKL  *  LKUX  r*Knr**2.0  *  6FR.ICK.)/ 
(66  *  DKL  *  tfHoKX  *  A*U  *  -*  2.0)  +• 
(6  *  £ke*  ^Hokk*  fc££A\x)  / 66 
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The  friction  factors  are  in  the  same  form  as  those 
for  the  downcomers 

FOE-  1.0/(l.74  “2  4>  &  ALOSI0  (ksDE)) 

FEF  -  1.0/ (l.?4  -Z0  #  AUP6I0  (fc^EF))  £rc] 

FJX  -  ±.0/(l.H4  -2.0  *  Ald6I0(ksTK^)  \Tca 

FKL  -  1.0/ (1.74  -2.0  x  ALOSI0  (ksRl))  \rcQ 


At  a  specified  steady  state  operating  point  the 
pressure  drop  across  the  downcomers  must  equal  the  pressure 
drop  across  the  risers  and  the  downcomer  flow  rate  must  equal 
the  riser  flow  rate.  Therefore  by  equating  the  appropriate 
pressure  drop  equations  and  solving  the  resultant  relation  for  the 
flow  rates,  the  initial  flow  rates  may  be  computed  as 

NF F00  -  ((KH0CD0  *  €v  -<a  *  XDE0  %.  ((ftV*QDO0 

-Fftl*OEE0)/2.0)—  6  *  £f=F0  *  RWOEE0*: 
<6iRA\/£)/((FCOjK  LCD/OCD  -FE>omCD  A* 

GEiODCO  -F  exiTCO)/(z.0^AC£)^^  2.0 
^  K.HO6O0)  -F-(,(rtHoe£0  -RHOOO0')/ 

(RWcttcfi  a  f?HcOO0  *  AOF  *  5k  z.0  ))  +- 
(4-0  5*  FDE  *  LDE0  *  7.0)/ (,7.0  *  DDF 

*  (*HOEE0  +^HODO0)-tK  AOF  *  *  2,0)  +* 
kACLE/C«Ho£E0  AOF^k  5K  2.0)  +  (^4-0 

*  FEF  5k  LEF0  #  £FK| CE.)/fe.0  3k  DOF 

*  R4OEE0  *  ADF  *  *  2.0)))  5k  *  0.5“  jrc] 
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NVU.00  = 

4-RHOKk0/2.^)-G,  H  ftKL0tt*H-OfcK0-* 
RSRAa/K*')  /  ((  F6H  *.  L<al4/06tt4  4-  EkJT^SH  +■ 
6EKJ06H  4-  E*IT6tt)/fr.0-fc  A^H  *  *  2-0 
^RttOS4^  -V-^KOKK^  -RHO3X0)/ 

*  RHOCTT0  *  ATL  *  *  Z.0))  +■ 

(4-9$  *  FT*  *  UTK0  *  Z.0)/(Z0  *  OTL 

*  (*HOkK0  4-  RHOJT^)  *  A<0-  *  *  2.$z0  +* 
RACxX/  (R.HOKK^  *  ATL  *  *  7.0)  + (q.0 

*  FKL  *  LKL^  *  RF=Ric K)/(z.^5  ^OTK 

**Hokk#S*  *tl**  *.£)))  **0.s-  \k3 


4 .  Riser  Continuity 

The  relationship  between  the  riser  inlet  and  outlet 
mass  flow  rates  is  written  in  terms  of  the  continuity 
equation  for  one  dimensional  unsteady  flow, 


In  model  notation,  this  becomes 

tALL  -  rv\TX-  draotl  *  voltl. 
nFF  =•  rADD-  ORH-OOP  ^  VOUDP 


Ime] 


The  numerical  differentiation  technique  used  by  CSMP- 
III  is  highly  suspect,  as  are  other  techniques.  An  "averaging” 
system  is  used  in  the  actual  dynamic  model.  In  addition,  the 
flow  rate  down  the  downcomer  and  into  the  riser  is  held  con¬ 
stant  for  the  open  loop  boiler  model. 
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5 .  Riser  Quality  and  Density 

The  average  density  in  the  risers  must  be  solved  for 
explicitly.  Linearly  varying  quality  along  the  riser  tube 
length  follows  directly  from  the  assumption  of  uniform  heat 
flux  along  the  riser  tube  length.  The  average  density  along 
the  tube  length  is  the  sum  of  the  total  change  in  density  and 
the  density  entering  divided  by  the  riser  tube  length  [1] . 
Assuming  the  density  varies  only  in  the  boiling  length  of  the 
riser,  the  total  change  in  density  can  be  written; 


where  =  density  of  saturated  water 

L  =  total  tube  length 


Since  p(&)  varies  linearly  from  entering  to  exiting,  p (  ?. )  can 
be  written  in  the  form 


\ /_£.  + 


*o3  (j-  lnj g  Vf*. 
L.s  T 


where  =  specific  volume  entering 

X  =  quality  at  riser  outlet 
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Lg  =  boiling  length 


! 

r 


LNg  =  nonboiling  length 

V  =  change  in  specific  volume  from  saturated 
®  liquid  to  saturated  vapor 

The  integral  of  the  explicit  equation  for  o(£)  may 

be  solved  after  rearrangement  to 

_ A- _ 

(v/jL  -  Ln/6  v/ly")  +  V-U 

Le  '  lq 

which  is  of  the  form  | 

a  + 

This  yields  the  average  density. 


f  In  1 

J-  +  1 

LM 

4- 


ft  L^_ 


In  model  notation  the  formula  is  written 


*HCdF0  -  (lEF^  /  ((XFF0)  # 

*  ALO€,(((XFF0)/\/F)^VF6i  4-1.0} 
4*  fcUO  000  *  LOE0  ) 


Ktt cxj\.<£  *  (i.0/UL^  (UO.^/CCXUL^')^  VF6?)^ 

Ac\_o6v  ((xui-^V/F6iyVF4'1.0)4-^)4c3^! 

)fc.  UK-0  £lC/A&] 
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To  permit  calculation  of  the  quality  term  the  energy 
balance  on  the  liquid  in  the  riser  tubes  is  evaluated.  Thus, 


rate  of  change  =  flow  energy  -  flow  energy 
of  riser  energy  in  out 

+  thermal  energy  in 

=  ^.vAv,  4-  ‘j.V. 

due 

If  the  enthalpy  term  is  treated  as  the  average  enthalpy  in 
the  riser  the  equation  may  be  written 

d  (s  (hf  +  x^*^)V)= 

In  model  notation  this  is  written  in  the  form: 

*  X  U-  *  FF6  /*•  <t>  )  (  (fcUOT L  *  VOlKL.”)  Od 

yUll  »  iure,Ru  (hxl\^}OHXi-i-)  E^e] 

^  (mu,  -  hf )^'Z.^/HF6» 

OvVXFF  «  *  (HOO-HF-XFF^ 

-r^FF  *  XFF  #  HF6/^)/(fcrtODF  *\/Ol£f)  [*Ej 
H*FF  =,  liOT^euCFXFF^jOBXFF)  [ V ie] 

XFF  -  (v+XFF  -  rtF)^^.^HF6i  tyc] 

The  initial  conditions,  HXLLO  and  HXFFO  are  calculated  in  the 
initial  condition  program. 

The  steady  state  equations  for  quality  are  derived 
differently.  The  energy  transfer  rate  to  the  riser  liquid 
may  be  written  in  the  form 

Q  =  rv\  A ,K 


where 


Ak 


=  V\  oui" 


—  h.fl  +  Xo«i* 


Solving  for  X  _  yields 
&  out  J 

*»>•*■  -  (§  Vw  (hf  -  li;*  ))  / (rA  kfy) 


or 

XFF0  «  (ql^  f^FF^  *  (HF  -HOO^)/(MFFyi  *  HF6.)  [icl 
X  Ll^  =  ($5#  F  Mil$S  *.  (i+F  -  HJ3V)  ft  /  (rtLLjzS  *  HFS)  [re] 


6 .  Riser  Boiling  Boundary  Location 

The  nonboiling  length  of  riser  tube  is  the  product 
of  the  riser  tube  length  and  sensible  to  total  heat  ratio  [l] 

The  sensible  to  total  heat  ratio  may  be  expressed  as  the 
enthalpy  change  in  raising  the  water  to  saturated  conditions 
divided  by  the  total  enthalpy  change,  i.e. 

9s .  -  k  -  k'n _ 

(hf  +  k^.)-  kv, 

The  equations  for  nonboiling  length  are 

Loe^  -  LDF  AmAXl^HF-vtoo / ((\\Qz>(f> 

Y  XFF£$  RFS)  -  MDO^  ) 

LTH^  IXL*  AmA*1/(HF-R 

-FXul^  *  HF$)  -  HCT3#j  tXlQ 


45 


The  Fortran  function  AMAX1  is  used  here  because  the 


above  equations  are  used  in  an  iterative  loop  in  the  initial 
condition  program  and  there  is  a  possibility  during  iterations 
of  reaching  a  situation  where  the  enthalpy  entering  is  greater 
than  saturation  enthalpy.  The  final  initial  condition  solu¬ 
tion  prevents  this.  The  AMAX1  function  is  not  used  for  the 
dynamic  model.  Thus, 

L0£  =  LOF  *  (hF  -HD0)/((  HF  -FXFF  #  HF<5i)-  HDD)  [we] 
UX  =  LTL  *  (HF  -  HTT)/(C  HF  *  HF^)  -HTT)  [wb] 


The  boiling  volume  is  derived  directly  from  the  solution  for 
the  nonboiling  length. 

boiling  volume  =  total  volume  *  boiling,  length  Hence, 

total  length 

VO l_EF  -  V&LDP  ^  UEF/LOF  |j*e] 

VOlKL.  -  VOLT*L  lkl/ltu 

7 .  Steam  Drum  Liquid  Mass  and  Energy  Balance 

The  rate  of  change  of  liquid  mass  within  the  drum  is 
equal  to  the  sum  of  the  mass  flow  rates  entering  and  leaving 
the  drum.  The  liquid  mass  flowing  into  the  drum  is  the 
saturated  liquid  from  the  risers,  the  liquid  from  steam  con¬ 
densation  in  the  drum,  and  the  incoming  feedwater.  The  liquid 
leaving  is  that  leaving  via  the  downcomers.  The  dynamic  model 
equation  is 


DOtfftSu  -  W»-L  *  (l.^-XuJ)  -hr^FF  *  (i.^-XFp) 

+  rv^CoiOJD  -H  mee  —  MCC  -  M<S  s 

The  instantaneous  mass  of  liquid  in  the  drum  is 

Om*SL  ~  l»JT6^U  (om*Su^  D-Dr^su)  I^e[] 

where 

DmftSL^  =>  (v&LDfc/m  yk  ^00e^u)/2.^  [xc] 

The  drum  energy  balance  is  derived  similarly,  i.e., 

DDWOHL  -  P\LU  *  (l 4  “XU.)#  HP  +  mFF  *  (l.^  -  XFp) 

*  +  mdDtJO  *  HF€»  +W3G  *  H6fc  - 

*  wu:  -  was  *  hss 


The  instantaneous  drum  energy  is 

Dm  DHL  -  l»JT<$^u(Drv\omj^>  Di^OH-l) 


The  initial  drum  energy  is  the  product  of  the  initial 
drum  liquid  mass  and  initial  drum  enthalpy. 


OmohljS  =  DmASl_0  *  HO^urvifS  \j£] 

The  initial  drum  liquid  enthalpy  is 

DH  -  Otad\\l i_ 

In  the  dynamic  model,  drum  liquid  enthalpy  is  considered  the 
enthalpy  of  saturated  liquid. 
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The  steam  condensation  rate  equation  is  based  on  the 
difference  between  the  pressure  and  temperature  of  the  steam 


and  that  of  the  liquid  [12] 

taCotJo  -  S<o<j>.  ^3  ^  (frM[  /(mtA  \L\(of.(j> )  tK  % 

-  fW  /  ( ■ TS  Arr  +  H  <4 .  <f> )  *  *.  .  r)  +-  $  2  S&  S' 


The  water  level  in  a  marine  boiler  is  generally  con¬ 
sidered  in  reference  to  a  mid-drum  zero,  that  is,  a  water 
level  of  plus  one  inch  implies  the  water  level  is  one  inch 
above  the  drum  centerline.  For  computational  efficiency  the 
assumption  is  made  that  for  small  changes  in  water  level 
around  the  midpoint  of  the  drum  the  water  surface  area  remains 
constant.  This  allows  a  simplified  level  equation,  i.e., 

LEv/Eu  -  (orAOV -  voud&*\  j% .<£  )  /  (Lsmoyv\  *  ostmjdm)  [me] 

DMOV  is  the  equation  for  the  total  volume  of  "liquid”  in  the 
drum.  Recalling  that  a  small  percentage  of  steam  leaving  the 
separators  is  entrained  in  the  liquid,  the  rate  equation  for 
"liquid"  volume  in  the  drum  is 

ODrooV  -  ((MFF  +  rvMJ--W\(3e)  %  VF  +PCi*  %  (tmll  -h^iff) 

w  ■+  >^66^  V66  — (fACd,  ]/F 

*i-rv\CDK)o  ^  l/F<s)  ]V\e] 


f 


where  DMDVO  is  the  volume  occupied  at  time  zero  which  is  half 
the  drum  volume. 

8 .  Circulation  System  Initial  Condition  Iteration 
Procedure 

The  initial  conditions  for  the  circulation  system 
are  found  by  flow  rate  balancing  in  the  downcomer-riser  flow 
loops.  An  initial  guess  of  flow  rate  is  made  using  the  assumed 
riser  exit  quality  of  .05  percent  which  is  reasonable  for  this 
type  of  boiler  [7] .  Coupling  this  assumption  with  the  assump¬ 
tion  that  the  percent  carryunder  is  zero,  an  initial  approxi¬ 
mation  of  the  flow  rate  can  be  determined  using  the  figures 
for  initial  heat  transfer  rate  to  the  risers  calculated  in 
Section  III-A.  Hence, 


mass  flow  rate  =  heat  transferee 


assumed  quality  *  latent  heat  of 

vaporization 


A  first  approximation  of  downcomer  enthalpy  for  both 
banks  can  be  calculated.  This  downcomer  enthalpy  is  assumed 
to  be  the  same  for  all  downcomers.  Using  an  energy  balance 
on  the  liquid  in  the  steam  drum 


MFP^  =  $^/(xASumE.  *hf€,)  Lrcl 

rALlfi  -  C^A-SuaiE  ri-Rs)  tX<£J 

^  + 
ttsfyrf  *  LxcJ 


The  screen  riser  inlet  enthalpy  is  assumed  to  be  equal  to 
that  of  the  screen  downcomer;  however  the  main  bank  fluid 
absorbs  additional  energy  from  the  desuperheater  located  in 
that  circuit.  Therefore,  the  main  bank  riser  inlet  enthalpy 
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must  be  calculated  separately.  Thus, 

4-  Lie] 

The  downcomer  density  must  be  calculated  for  use  in  the  pres¬ 
sure  drop  calculations.  Hence, 

VCo^  -  ((rwFF^  \A=  4-  / 


tic] 

RH Qtctf 

XF<Q 

- 

The  riser  outlet  quality  is  calculated  using  the 
initial  quality  formulation  previously  developed.  This  com¬ 
putation  is  followed  by  the  calculation  of  the  riser  non¬ 
boiling  and  boiling  lengths, 

The  average  density  of  the  risers  is  calculated  along 
with  the  two  phase  flow  multiplication  factors  for  use  in  the 
flow  rate/pressure  drop  calculation. 

An  updated  flow  rate  is  now  computed  and  compared 
with  the  first  approximation.  If  it  is  within  a  specified 
error  criteria  calculation  stops.  If  not  the  previous 
approximation  is  updated  and  the  calculations  continued  with 
the  new  approximation. 

Upon  completion  of  the  flow  rate  balancing  the  initial 
flow  rates  in  the  downcomer/riser  loops  are  known  along  with 
the  riser  outlet  quality. 

'the  steady  state  drum  specific  volume  and  density 
are  computed  as 
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\JOR. -  \f( 100  4-  pOu  *  \}\j 

RoofcmL  -  l.p 


The  initial  steam  mass  in  the  drum  is  then  calculated  as 

0STfv\<jz$  -  voi_o e^\  4K  £rto\/ 

9 .  Superheater  Pressure  Drop 

All  entrance,  head,  and  exit  losses  are  considered 
negligible  in  the  superheater  compared  to  the  frictional 
pressure  drop. 

iP  =  |  L.  i^Z 

A2-  ^  averse- 

By  lumping  the  constants  the  pressure  drop  equation  is 

Pwv\-(°M|0  ~  tCOKJl  tK/VW*!  ^  ^  ^  ^  j 

The  following  definitions  apply  here: 

Ko0  3-  (pr*M0  +  pMD^)/(f^H-onrn^-VRVIoN3NJ0)  ^3 
KHomnJ^  -  (RHowrvi^  4- 

KoiOA.  —  ^ 


10 .  Steam  Valve  Equation 

The  flow  through  the  steam  valve  is  considered 
directly  proportional  to  the  outlet  pressure  and  valve 
opening,  i.e., 
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M  —  C/  ^  P  ( v/aWe.  Ope^ina^) 

Mwnz  -  pmiO  *  KofJ4  *  valv/e.  \Vv\e] 


The  constant  is  computed  in  the  initial  condition  program. 

Ko\)4  -  irwrwmxxi:  / (vau/e$  *.  PnJnJ^)  \jccl] 


E.  EQUATIONS  OF  STATE 

The  equations  of  state  listed  below  are  used  in  both  the 
initial  condition  and  dynamic  programs.  With  the  exception 
of  the  subcooled  specific  volume  equation  used  for  the 
feedwater  entering  the  drum,  they  are  reasonably  accurate 
in  the  300-1500  psi  range.  The  subcooled  specific  volume 
equation  is  accurate  in  the  £oo  -  S'oo  psc  range . 

PSAT  -  EXP((AU>$  4(*7^*)/.2<W^  ]X 

TSAT  EXp((  .22.15*1  ^  AbCA(p5AT)  +4.H7lVo))  K 
HSAT  -  EXP  (0.26452,  *  ALOS  (psat)  +•  4.4670?*)  M 

ft/AT-hi.7i7e-^4 

&PSAT*  *  2.0  —  4.2.14E-0?**  3-0  I<1 

RHCF  -  (dZ>%  -0.0i7?l  #  TSAT -V  1. 1*32.  E. -05* 

*  TSAT  *  *  2.0  -  *  TSAT 

*  *  3-0  6] 

V66  -  .0l(o004?S  -  •  0^00^2.0  H(c>  #TB6 

+*  3.  Grua  -0?  *  tss  *  *  2.0 
-  ?.142E-II  *  TSG  F1.4^\e-i3 

*  Tfcfl  *  4.0  -  UI4?E-\(p  *-03(2**  S*.0 
■H  ?.C54E-20  *T3G  **6.0 

i 
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IV.  RESULTS  AND  CONCLUSIONS 


A .  GENERAL 

A  listing  of  the  initial  condition  program  and  the  dynamic 
boiler  model  program  is  given  in  Appendices  A  and  B  respec¬ 
tively.  The  initial  condition  program  output  data  must  be 
properly  formatted  for  input  to  the  CSMP-III  dynamic  model. 

In  addition,  because  the  dynamic  model  utilizes  the  liquid  in 
the  steam  drum  as  a  saturation  state  point  from  which  all 
other  state  points  are  derived,  two  of  the  initial  conditions 
must  be  modified  to  eliminate  a  discontinuity  between  the 
steady  state  program  and  the  dynamic  model.  The  equations 
involved  and  an  explanation  are  given  below. 


ompiSl^  -s  (vold*w  *  Rt+c&OfO  Irc3 

Dn'iO  H  \-<f>  =  »4c*-um0 

£>/vU>H-L  ^ 

DlY)f\ SL  -  \t)rG&L.(DtOfiSL0  LmS0 

OVV  -  OmDHL./o^AS’L 


As  stated  previously,  the  enthalpy  of  the  liquid  in  the 
drum  is  considered  saturation  enthalpy.  The  dynamic  model 
must  begin  with  the  rate  of  condensation  in  the  drum  (MCOND) 
as  close  to  zero  as  possible.  This  is  a  natural  steady  state 
pos it  ion  -  drum  liquid  and  drum  steam  both  at  the  same  pressure 
and  temperature.  In  order  to  insure  this  to  be  the  case 
DMASLO  and/or  DMDHLO  must  be  modified  such  that  the  initial 
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drum  enthalpy  (DH)  very  closely  approximates  the  enthalpy  of 
saturated  liquid  corresponding  to  the  drum  steam  temperature 
and  pressure  (PMM  and  TMM).  This  is  easily  facilitated  by 
the  use  of  the  CALL  DEBUG  statement  in  the  dynamic  program 
With  this  procedure: 

1.  The  initial  condition  program  is  executed  and  the 
initial  conditions  formatted  for  CSMP-III  use. 

2.  The  dynamic  model  is  executed  for  only  a  short  run 
time,  i.e.  5  seconds. 

3.  Observing  the  DEBUG  output  from  the  model,  DMASLO 
and/or  DMDHLO  are  modified  such  that  PSAT  and  TSAT 
equal  PMM  and  TMM. 

4.  The  model  is  reexecuted  for  a  short  run  to  check, 
MCOND  should  be  very  small. 

5.  Since  DH  is  implicitly  related  to  PMM  and  TMM,  the 
procedure  may  have  to  be  repeated  a  few  times.  The 
objective  is  to  force  as  many  of  the  "DERIVATIVE" 
terms  in  the  DEBUG  output  to  comparatively  small 
figures  as  possible. 

The  model  as  used  operates  satisfactorily  with  DMASLO  = 
8002.2  and  DMDHLO  =  3.96163E  06.  This  forces  the  initial 
MCOND  term  to  .35949. 

Because  of  the  CSMP  function  DERIV  used  in  the  program, 
the  model  is  extremely  sensitive  to  integration  time  step. 
Numerical  differentiation  is  not  a  desirable  function  to 
perform  in  a  dynamic  model;  however  attempts  to  explicitly 
differentiate  the  equations  concerned  failed.  One  solution 
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to  the  problem  is  to  "smooth”  the  derivative  function  by 
averaging  it  over  several  timesteps.  The  model  performed 
satisfactorily  with  a  fixed-step  integration  procedure 
( Runge-Kutta) ,  an  integration  time  step  of  .04  seconds, 
and  averaging  the  derivative  over  sixteen  time  steps. 

B.  OPEN  LOOP  RESPONSE 

The  open  loop  response  to  a  ten  percent  increase  in 
throttle  opening  is  sho-vn  in  Figures  2-7.  As  expected,  the 
response  of  the  main  bank  circulation  loop  and  screen  bank 
circulation  loop  is  different.  This  is  in  keeping  with  the 
■  different  purposes  of  those  two  loops,  steam  generation  and 

furnace  screening  respectively.  The  effect  of  the  valve 
opening  increase  is  barely  noticeable  in  the  screen  circula¬ 
rs*  tion  loop. 

The  swell  effect  is  not  noticeable.  Further  conversa¬ 
tions  with  Mr.  Paul  Weitzel  at  Babcock  and  Wilcox  indicate 
that  the  percentage  of  "carrv  under"  is  not  a  constant  as  it 
was  treated  in  this  program.  One  to  two  percent  is  a  good 
starting  estimate  for  steady  state;  however  during  a  transient 
,  "carry  under"  mass  flow  rate  is  computed  by  subtracting  the  amount 

of  steam  leaving  the  boiler  from  the  amount  of  steam  produced. 

This  should  always  be  a  positive  number.  The  vapor  that  does 
not  leave  the  boiler  is  "carried  under."  The  following  program 
will  implement  this  "carry  under": 
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Modify 

the  section  of  the  dynamic 
model  titled  "COMPUTE  THE 
DRUM  SPECIFIC  VOLUME"  to 
read: 

PROCEDURE  CRYUND=F I LTR8 ( MMM ) 

IF(TIME . GT . 0 . 0 )G0  TO  45 
CRYUND=PCU* ( MLL+MFF ) 

GO  TO  46 

45  CRYUND= (MLL*XLL+MFF*XFF ) -MMM 

46  RISE1=DELAY  ( 250 , RISTIM , CRYUND ) 

RISE2=CRYUND 

IF ( VALUE. GT. . 51)G0  TO  55 

RISE=RISE2 

GO  TO  57 

55  RISE=RISE1 

57  CONTINUE 

ENDPROCEDURE 

DDMDV= ( ( MFF+MLL-MBB ) * VF+CRYUND* VV . . . 

+MBB0  * VBB- ( MCC+MGG ) *VF+MC0ND . . . 

*VFG-RISE*VV) 

Follow  the  above  statements  with  the  unmodified 
equations  for  DMDV. 

The  previous  procedure  was  not  implemented  in  the  present 
dynamic  model.  The  quality  formulation  in  the  present  model 
is  apparently  too  simplified  and  this  causes  a  shortfall  in 
quality  at  the  outlet  of  the  screen  riser  bank.  As  a  result  of 
this  shortfall,  there  are  instances  during  transient  operation 
when  the  model  is  not  "producing"  as  much  steam  as  it  is  "using". 
Mr.  Weitzel  suggests  using  twenty  node  finite  difference 
approximation  for  quality.  This  could  be  done  using  an  equation 
for  quality  such  as: 

x  -  ...  q _  * 

A  Si 

where 

X  =  quality 

Z  =  distance  up  the  riser  (ft) 
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=  heat  input  to  riser  (BTl'/s-ft) 

2 

A  =  cross  sectional  area  of  riser  (ft  ) 

3 

c j  =  density  of  saturated  liquid  (lbm/ft  ) 

=  velocity  of  liquid  entering  (ft/s) 
h^.  =  latent  heat  of  vaporization  (BTU/lbm) 

Simultaneously  an  implicit  equation  for  continuity  could  be 
applied  to  produce  a  balanced  mass  flow  rate  without  assumin 
that  the  downcomer  mass  flow  rate  remains  constant. 

C.  CONCLUSIONS 

The  model  presented  is  not  a  finished  model.  Further 
research  is  needed  to  successfully  implement  the  shrink  and 
swell  theories  presented.  When  the  previously  mentioned 
difficulty  with  riser  outlet  quality  is  solved,  the  model 
should  approximate  the  dynamics  of  a  wide  range  of  D-type 
marine  boilers  depending  on  the  initial  conditions  supplied 
to  the  initial  condition  program. 

The  initial  condition  program  develops  the  initial  condi 
tions  necessary  for  the  detailed  boiler  model  with  the 
relatively  scant  data  available  from  the  boiler  technical 
manual  shown  in  Appendix  C  and  a  small  amount  of  data  from 
common  engineering  handbooks.  Because  of  the  relatively 
small  amount  of  data  required  the  initial  condition  program 
and/or  the  model  can  serve  as  a  basis  for  a  boiler  condition 
monitoring  system. 


D.  SUGGESTIONS  FOR  FURTHER  RESEARCH 

It  is  evident  that  the  model  suffers  by  using  the  CSMP- 
III  function  DERIV.  In  order  to  facilitate  the  explicit 
solution  of  the  continuity  equation  for  the  riser  banks  some 
form  of  differentiation  is  required,  either  numerical  or 
explicit.  An  explicit  solution  is  much  more  desirable  from 
a  stability  standpoint.  A  more  detailed  investigation  with 
fewer  assumptions  might  produce  the  correct  form  for  the 
explicit  differentiation  equation.  It  should  be  noted, 
however,  that  when  using  explicit  differentiation  the  model 
must  average  the  derivative  over  at  least  two  time  steps  to 
avoid  the  creation  of  an  algebraic  loop.  This  is  easily 
done  with  the  same  PROCEDURE  format  used  to  average  the 
derivative  in  the  current  model. 

The  use  of  small  perturbation  techniques  to  linearize  the 
model  and  thus  allow  a  state-space  representation  for  multi- 
variable  control  development  and  analysis  should  be  undertaken. 
There  is  a  computer  code  available  locally  to  facilitate  this 
for  CSMP  models;  however  it  does  not  accept  the  CSMP  DERIV 
function  and/or  the  averaging  procedure.  A  means  of  bypassing 
this  problem  would  result  in  a  general  D-type  boiler  model  - 
applicable  to  a  very  wide  range  of  boilers  currently  in  use  - 
which  could  be  linearized  about  specific  operating  points. 

These  individual  models  could  then  be  used  for  the  development 
and  testing  of  multivariable  control  systems. 

Using  locally  available  optimizing  computer  codes,  an 
optimal  D-type  boiler  design  could  be  attempted  with  regards 
to  boiler  geometry. 
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A  more  detailed  analytical  and  experimental  investigation 
of  current  D-type  marine  boilers  should  be  undertaken  locally. 
Little  research  has  been  done  in  this  area  and  that  which  has 
been  done  has  often  been  based  on  either  an  incorrect  physical 
model  of  a  marine  boiler  or  on  a  nuclear  steam  generating 
plant.  A  starting  point  could  be  the  data  analysis  of  common 
horizontal  and  vertical  steam  separators  followed  by  an  opti¬ 
mum  design  for  these  elements.  The  improper  and/or  mainten¬ 
ance  of  drum  internals  apparently  grossly  affects  boiler 
operation  and  water  level  stability  during  rapid  transients. 


TABLE  1 

MODEL  NOTATION 


NOTE:  In  general,  rate  variables  gr  are  preceded  by  the  letter  D,  i.e.,  DHXFF,  DRHOJL. 
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£8/3075  2C.55.II 


FILt:  CONSTANT  FORTRAN 


NAVAL  POSTGR AQLATE  SCHCCl 


c 


this  progra-  calcllates  ’■he  constants  anc  initial  cdncitidns 

FOP  US  NAVY  D-7YPE  H  JI  Lr  R  GIVEN  THE  INPUT  JATA  OBTAINED  fRC" 

THE  ECILES  TECHNICAL  MANUAL,  STF.V*  TABLES,  AN 0  STArDLR)  ENGINEER¬ 
ING  HANOBCUnS.  THE  INPUT  CAT  A  ARE  ENr  ER  EO  IN  N A  A  E  LIST  FORMAT. 
NAMELIST  DESIGNATIONS  AND  INPUT  NEWONICS  FJLLDw 


INCQNl: 

OPPNT  OPERATI  JG  PO  I  NT  I  PERCENT )  * 

TOTSTM  TTTAL  STER*  FL3*lLB/HR) 

SFhSTN  SUPERHEATER  STEAV  FlCW  (LB/HR) 

0SH3TN  DF3UPERHC ATEO  S  T£  A  M  FLOW  (LB/HR) 

PCRUM  CSLM  PSESSURt  (PSIG) 

SHOT  SUPc  Rh  EA~£  R  CUTLET  TEMPERATURE  (CEG-F) 

SHOP  SUPER  FEATt3  OUTLET  PRESSURE  1?SIG) 

OSHCT  OESU  p  EPHE  AT  ER  IUTLcT  T  £mp  E°  ST' jo  E  ( JEG-F  ) 

OSHDP  DESUPERHEATER  OUTLET  PRESSURE  IPSIG1 
ECONIT  ECONO'*  12  E3  rEE  1  INLET  TEMPERATURE  (JEG-F) 

E CONST  ECCNCM2EF  FcEJ  OUTLET  t  Ey  oE  *  yr  :J 0  E  (CEG-F) 

A  I®  T  HP  AIR  INLET  TEMPERATURE  TO  REGISTcRS  (JEG-F) 

CILTVF  CIL  T  EvPfcRATURt  AT  =J3NER  INLET  I JE G-F  ) 

AIRFLO  AIRFlOW  RA~E  (LB/HR) 

0 ILFLC  OIL  FlO*  RATE  (LP/HR) 

XCSAIS  £<  CtS  S  AIR  (PcRCE'D  * 

DP  AFT  DRAFT  LJSS  (  IM-H2Q) 

TC-A3SC  FLUE  GAS  TEMPER  Mu®  F  LEAVING  SCREEN  (DEG-F) 

TGASSH  FLLu  GAS  -EMPtRATjPE  LEAVING  SU°E3HEATE3  (CEG-F) 

tgasmp  flue  gas  temperature  leaving  -ain  bank  (jeg-fi 

TG  AS  EC  FLUt  l  AS  temperature  LEAVING  £r  JN  )M  1 2  ER  (DEG-F) 
TSCRN  SC  RE  E\  TUBE  WALL  TEMPERATURE  (DEm-F) 

TSPHTR  S;J°  £R  F  EATER  TUBE  WML  TEMPERATURE  (JcG-F) 

HkFVCL  Ht  Ar  R  cl  £  i  S  £  (<?TU/fr/cj  ft  F'JONACS  VOLUME) 

FHAS  FURNACE  HEAT  ABSORPTION  (KBTJ/HR/SQ  FT  R  A  C  I  ANT  HEAT 

ABSORBING  AREA  ) 


NOTE.*  "PERCENT"  VALUES  ARE  TO  BE  LCGGED  AS  DECIMALS. 

*  SCREEN 

DTUBSC  A'/ERACE  1  *S  IDE  DIAMETER  OF  Sro  EEN  TUBES  (IN) 

LAV SC  AVERAGE  LENGTH  GF  SCREEN  TUBE  (FT) 

NTUBSC  NUMBER  IF  SCREEN  TUBES 

RHASSC  SAJI  Ai'.T  HEAT  ABSORBING  AREA  HF  FJRNACE  SCREEN  (SO  FT) 
MASSSC  TOTAL  WEIGHT  OF  SCREEN  TUBES  I  LB) 


JSPHTR 

AREASH  TOTAL  SUPERHEAT  ER  HEAT  TRANSFER  AREA  (SO  FT  ) 
MASSSh  TOTAL  WEIGHT  OF  SUPERHEATER  TUbES 


$MNBANK 


OTUBM  B  AVERAGE  ISSIOE  DIAMETRo  r,F  main  34,JK  TUBES  (IN) 
LAVM9  AVERAGE  L  ENG’ H  FF  MAIN  BANK  TUBE  (IN) 

NTU3M 9  NUM3ER  OG  MAIN  BANK  TUBES 

►ASS"B  T  CT  AL  WEICfiT  nG  m  4  Ra-lK  TU°ES  (L?) 

AREAM9  TOTAL  HEA T  TRANSFER  AREA  (SO  FT) 


1ECON 

DTLBEC  .INSIOE  0  I  A  VE1"  ER  CF  ECONOMI2ER  TUBES  (IN) 
NPASSE  NUMBER  0 F  TUBE  PASSES 
NTU9  EC  MuocS  OF  TUBES  0  E3  PASS 
L TU9EC  AVERAGE  LENGTH  OF  ECONOMIZER  TUBE  (FT) 
MASSEC  TMAL  WEIGHT  CF  EOONTMIZER  IlB) 


CONOOO  10 
CCA  TOO  20 
CONCCCA J 

CDNGuuhJ 

C0N00050 
CON  CCCfcO 
CON  0  JO  70 
CONOOO  HO 
CON  C  CO  SO 
CONOOIOJ 

coNOono 

CON  00  I  20 
C0N001  30 
CONCCI AO 
CONCC  150 
C3NC316J 
C0NCC1 7C 
coNOoiao 
C3NC01G0 
C0NC02CC 
C0N0J2 10 
CONCC22C 
C3NC02JJ 
CCKCC2 ,J 
C0MCC220 
C3N0J2o0 
C3N00270 
C3NCG230 
C3N0029J 
C0NCC3C0 
CONGO  3  10 
COMOO 020 
CDNCC330 
CO  N  0  0  G  AO 
CONC035O 
COf)CC3tO 
C0N00370 
C0NCC3  EC 
CON 00  3  AO 
C0N004G0 
COM  CCA  10 
CON C 04 2 J 
C0N30420 
CDNCC44C 
CQN0G457 
CONC0440 
C0NCC470 
C3N00460 
C0NCJ49D 
CON  00  500 
CO 4 C 05 10 
CON  CC  520 
C0N00530 
CON00540 
C0NCC550 
C  0  NOO  5  oO 
C0MCC570 
CONGO  5  FO 
CDN00590 
CON  CCtCC 
C0N0J6  10 
CDNCOfc  20 
CON  C  C6  30 
C0NJJ640 
C0NCC65J 
COM  C0660 
C0N006  70 
COM  CCfccO 

Con  ccg so 

CONGO  700 
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FILE:  CONSTANT  FORTRAN  PI 


NAVAL  POSTGRADUATE  SCHCOl 


$CES°H 

0TUB3S  INSIDE  HAMSTER  OF  DESU°£RhEA  TE  R  TUBE  UN) 

MUSCS  NUMBER  IF  UK  J  °  ER  HE  \T£R  TUQES  PER  J  ASS 
NPASSJ  NUMBER  CF  DESUPERHEATER’  TUBE  PASSES 
LTIIBCS  LENGTH  IF  0  E  SU°  5  R  HE  \  TF  5  TLhE  (FT) 

APEAC5  TTTAL  heat  TRANSFER  A3  E  4  OF  0  ESU°E  RH6  AT  E3  (SO  FT) 
MASSOS  TOTAL  WcIGHT  CF  De  SLPGRHE  ATER  ASS  6MB LY  (Lei 


1DRMDCR 

OTU30Q  AVERAGE  D IAMETE  R  CF  DRUM  CCWNCOMER  TUBES  (IN) 
LAVCC  AVERAGE  LENGTH  OF  OR  UN  DOWNCOMER  TUBE  (FT) 
NTLBCC  NUMBER  OF  CRUM  CCWNCOMER  TUBES 


CONCC71U 
COM007 20 
CD.N00  7  JO 
CON  C  C  7  40 
CONDO  7  5 C 

cancel  ho 

CUNCO  77G 
C0ND073U 
CONCC7SC 
CDN0030J 

concob  u 

CDNCC620 
C3N008  30 
C0NG0840 
CON  00  8  50 
CONCCB60 
C )N  CCe  7C 
CON JOB  30 
C  ONC  0  3‘>o 


SHDROCR 

0TU9HC 
LA  VHD 
NTU8H0 


JBCILER 

0  STM DM 
LSTMCM 

OWTRDM 

LWTRJM 

HN'CRM 

hhor 

h*T.RCM 

FLRVCL 


AVERAGE  INSIDE  DIAMETER  Hp  HEADER  DJWNCJMER  (IN) 
AVERAGE  LEAutH  CF  HEADER  C'NUCCIM  £R  (FT) 

NUMBER  JF  HEADER  DO  '«NC 0  ME  T  S 


DIAMETER  OF  STEAM  DRUM  UN) 

LENGTH  RF  S t E A M  DRUM  (FT) 

DIAMETER  CF  WATER  ORUM  (IN) 

LENGTH  JF  WATER  JRU«  (FT) 

He  1  vhT  :f  \ 0RM A L  V.5TET  LEVEL  ABOVE  BENCH  ma^k  (FT) 
HEIGHT  CF  HEADER  (SCREEN)  ABOVE  BENCh  map*  ( FT ) 
HEIGHT  OF  DRLm  ABOVE  BENCH  MARK.  (FT) 

FLRNACE  VOL  IME  (CU  FT  ) 


NCTE:  CHOICE  OF  BENCH  mark,. is  ARBITRARY 


(THERMO 

H  SHOUT 
HOSOUT 
HECIN 
HECUUT 
(1JKH20 

(1 )  FRH2  c 

(  l  )  VSC  H  20 
(2  >KSTm 

( 2 )  PR  STM 
(2  JVSCGTM 

R  SHCL1 
ROSDUT 
RFL'JE 

NOTES: 


SLCSSES 
KSCSC 
KSDMB 
KSCCC 
K5DH0 
ENT  SC 
EENCSC 
ENTMB 
B  £NCM  6 
E  MOD 
RENDOO 
ENTHC 


ECONOM  I  ZE  R  CUTLET  ENTHALPY  (  RTU/L  BM  ) 

DE  SUP t  RHE  A  TER  OUTLET  ENTHALPY  (3TU/L0M) 

ENT  h  A  L  PY  CF  6CONOMI265  FEED  INLET  (R7U/LB«) 

ENTHALPY  OF  ECONOMIZER  FEED  OUTLET  (ETU/LEM) 

therm ul  conductivity  j=  water  ibtu/hr.ft.deg-f) 

PRANDTL  NUMBER  FC  R  N AT  E  R 

KINEMATIC  VISCOSITY  FOR  WATER  (LBm/FT-SEC) 

^hERMAL  CONDUCTIVITY  FOR  steam  (  BTU/HR  .FT  .DEG-F  ) 

PRANDTL  NUMBER  FOR  STEAM 

KINEMATIC  VISCOSITY  FOR  STEAM  (LBM/FT-SEC) 

SUPERHEAT  ER  D'JT  LET  DENS  tTY  (LEM/C1  FT) 

DE5UPERHE ATER  ULTLET  DENSITY  (LBM/CU  HT ) 

CENSITY  of  FLUE  GAS  AT  TGASSC  (LBM/CU  FT) 

(1)  EVALUATED  67  AVERAGE  ECONOMIZER  WATER  TEMPERATURE 

(2)  EVALUATED  AT  AVERAGE  DESUPERHEATER  STEAM  TE  MPER  A  TUC 

ROUGHNESS  RATION  (S  AND  EQUIVALENT)  FOR  SCREEN  TUBES 
ROUGHNESS  RATIO  (SAND  E  OUI VALENT)  OF  MAIN  BANK  TUBE 
ROUGHNESS  RAT  U  (SAND  EQUIVALENT)  FDR  D R LM  DDWNCOwE° 
ROUGHNESS  RATIO  (  5- NO  EQUIVALENT!  F CR  HEACER  DOWNCOM£R 
SC  3  6  E  U  TUBE  ENTRANCE  LOSS 
S  CREEN  TUBE  0E‘<O  L~SS 

main  ba  jk  entrance  loss 

MAIN  BANK  BEND  LOSS  FACTOR 
DRUM  CC  VNCCMER  E  NT  R  ANCE  LOSS  FACTOR 
DRUM  J  CWNC  CiMc  R  BE  ND  LOSS  FACTOR 
HEACE.R  CCWNCOMER  ENTRANCE  LOSS  FACTOR 


CONOOSoO 
CCNJJ9 IJ 
CONG  OS  20 
C0UGGS3J 
CONG 0940 
CON  CCS  50 
CCN00960 
CDN0C97D 
CON  CGSRO 
CQN00SSO 
CONCIOOJ 
CON  CIO  10 
CQNJIJ2 J 
CONC  ICC 0 
CONG  10  40 
CONG l J  50 
CONCICfcO 
CD  ND 10  70 
CON 01 CBO 
CDNC1CS0 
CONOUJJ 
CONClllO 
CDNCII20 
CGN01I30 
CONC  l  140 
CDMCU5J 
CO NO  1 160 
CONC1170 
CONO 1 180 
CONO  1  190 
CONO 1200 
C0MC1210 
CON  0  12 iv 
CONO  1230 
CON  C  1 2  40 
COMO  12 50 
CONO 1260 
CONC 12  70 
CONO 12  30 
C0ND1290 
CONC  13C0 
C0NO1310 
C3NC1 320 
CONO  13  JO 
CO  NC 1340 
C3NC1350 
CD  N  0 1360 
CONC1370 
CDN0139G 
CONO 1390 
CONC  14C0 
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FILE:  CONST  ANT  FORTRAN  PI 


NAVAL  POSTGRADUATE  SCHCCL 


c 

c 
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8EN0H0  HEADER  DCMNCCMER  bend  l^ss  factor 

IMPLICi"  PE  AUA-2) 

INTEGER  NPASSE 

DIMENSION  Th23(  2G).QPASSI2C) 

DEFINE  NAMELISTS 

NAMEL I  ST  /  INCCNO/CFPNt  , "07S TM,S fhS t m,  CS  HS  T”  ,  PORU M,  S  HC T ,  SHOP  ,  0 SF^T 
IDS  FOP,  EC  Oil  T  » ECON  )  T ,  MR  T MP  ,  J I  LF  L3,  A  I  P.F  LC,  XC  S  41  R  ,  UR  A  F  T  ,T  GAS  SC  , 
2TCASSF,"  CAS  m*."  CASES  ,T3C  M,  TSPHTJ  ,  HP  FV"U,  Fha  0,  HL  IMP 
NAME  LI  ST/SCRt£N/CT'jaSC,LA7SC,NTUBSC,  SHASSC,M  ASSS  C 
NAMEL  IST/SPHTP/ARFASH,  -'ASSSH 

NAMEL  l"T/  v?;t  ANK./  LT  I  bvB,  L  AV  «3i'ITUS«a,  MASS’*  A,  ARE  AMR 

NA  ME  LI  ST  /ECCN/JTLJEC  ,NP4SS  E.NTUfcEC  ,LTH3fcC,  MASS  EC 

NA«  EL  IST/CES’  H/3TJ  VDSUTCPOS.’IP'-  SS3  ,  LTU3DS  ,  A  RE  AO  S ,  “  A  S  SD  S 

FAFELI  JT/Di,mDCR/LT-JPCC  ,CAV  CO,  MUR  03 

NAMEL  1ST  /r*3R  OCR  /  3  TiJ  UO  ,LA  VHJ  ,NTU3N J 

NAFELIST  /  5JIL  ER/  OSr-’iCM,  LST  ■■'j-*,  JwT1-'  O'* ,  L  WT R DM  ,  rNOR  M  ,  HHDR  , HWTRD M  , 
SFLRVCL 

NA«£L  IST/THE-'--n/HS  r>UT,HOSDUT,  HSCIN  .HECOUT, 
lKH2C,PRn2C,'/;  C.-i2  C  ,  <ST«,  PRST'1  ,V$rS7M,R  SH"'1  IT  ,  ROS  O’JT  ,  RFLUE 
NAMEL  I  sT/LOSSES/nSOSC  » Kc  0**6  ,is3JJJ,KSOtJ  ,EN7S  C  ,  &  EK  CS  C  ,  ENT  *3, 
ABENCM6,  E  ITCO,  3E  I  CC  J,  fN  THJ ,  0  END  HO 
NAMEliST  / 1 NCCNl  /”RRD  ,7  7V0  ,~V,WJ  , 

ETYYC,TZZ  C,  tXAJjTA'J 
NAMEL  IS  r  /  l  \  C  CN2 /  '•' A  VO  *  MFQQG  » M  AOOQ, 
fi  NNNO  ,MMMC , M JQO , 

EMCCO.MGHO.vPVO, 

E*FF 0 .MLLO.Wt, 

EMASSV*MAJS'«»  MASSY* 

£M  ASSX',  tASSZ,  -'ASSOR, 

COPUHIC, JHASLO, JSTVO  ,CSTMDM, 

EL57MDM , JST  10 M 

NAMELIST/ l  N CON 3/ rvN, CA3.CPV, 

ECPY  ,CPX  ,CPZ, 

ECOR.CQR  , CRS ,CST , 

£  CTU  ,C  N  P  » G »  GC 

NAME  LI  ST /CONST  1 /N SJL ,KSOF , K5GH , 

£KSCC,KSOE,KSEF, 

EKSJK  ,KSKL ,N2 , 

EKW  ,K°S,KJf, 

EKQM  ,  KO'I  j  , K  CN4 , 

£KX,KHJ,K\, 

EKT  U  ,  K  Y 

NAMEL1ST/CjNSt2/ AJL, CJL,LJL2, 

ELJLl.DJF ,L3F , 

SLCFl, ACF,L CF, 

EACH  ,Cl»H»CCJ  , 

ELC  D  ,LO£  0 ,  l  J  .<  C 

NAMEL IS" /CONST  3/ ACC, 23 ENCi.Z 3 cNC 2, 

E  VOL  J  L  ,V  C  l  )F  .VClOhM  , 

EVOL  H J , 2  OF , 2  JL , 

S LJL  ,CDF ,JKL , 

ECEF , 0 JK , 2CO 

NAMELIST/CUNST4/ FCO, FCH, FOE, 

EFEF  ,F  JK  ,F  A L  , 

♦ENTRGH.E  IYP  CO  *  3EN0GH, 

ebencco.exi'gh, ex itcc 

NAME  LI S  t/C  1N3T3/FHV  ,  XLLO.XFFO, 

CT  A*8 ,  C  5  0  CFG,  OR  JJLU, 

ESIGMAA 

NAMEL  l  ST /OUT3  L  T/C  10  ,C40 ,C 6C  ,OS C ,Q9 0 , MFF 0 , MLLO, CCO , 

EXFFO  ,XL  10 ,"  HCJFO  ,  CHOJLJ,  RH'OF  F  0 , 5  HOLLO, 

EHSAT  fHCOC  ,HAAO  ,H  8  30  ,hJJ0  ,HFG  ,  RhCF,  D*CVO  ,  HX  FFO ,  hXLLO 

/  Rc  AD  15  ,  INCCNO  ) 

READ! 5, SCREEN) 

RE  AC  15, S'1  FT3  ) 

•  READ (5, MN8ANIO 


C3N0  14  10 
C  ONO 1420 
CONC142C 
CON  C  1  4  *,0 
C0N0U5J 
CON  C146C 
CONOI4  70 
C  0  N  0 1 4  8  0 
CJ\C  US  J 
CCNOISOO 
C  ONO  1510 
C  ON  0  1 5  2  0 
CGN015 JO 
CON C  15  40 
C0NGI550 
C0N0156J 
C0NC157G 
TONCISBO 
C0NC15S0 

cor.'cicco 


C0NCI6 10 
COM C  1620 
COf'O  lb  30 
C  0  N  C  l  b  4  0 
CO  N  C  1  fc  5  0 
CONO 16  60 
C0N01A7C 
CON  016 80 
CCN01690 
C0NC17C0 
CON  C  1710 
CONO 17 20 
CONCl 730 
C0N0174O 
CONCl  7 50 
CONCl  7 tO 
CON01770 
CONCl 7EC 
CONJ 17SJ 
CONCl 8  00 

coNCiau 
C0N0U20 
CONOl 83C 
CON  C  18  40 
CONCl 850 
CONO  I860 
CON  0 18  70 
CONCl a  30 
C0NC16S0 
C0N01900 
CONO  19  10 
C0NC1S2J 
C0F0193J 
CON  019  40 
C0N0195J 
CON  C  1960 
CON C 19  70 
CONO 1980 
C0NC1SSC 
03N02000 
C0NC2010 
CON  C  20  20 
CONO 2 030 
CONO  2  C4C 
CON  C  20  50 
CCN02060 
C3NC2C70 
C0NC2C80 
CUN02U90 
CON  C2  1  CO 
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FILE:  CONSTA'r  Fortran  Pi 


NAVAL  POSTGRADUATE  SCHOOL 


c 

c 


REAOt  5,ECCN) 

SEAC 15,  C  tSn  H  ) 

READ (5,0* 'JCP) 

REACH  5,H  }A,)C  5  ) 

PE  AC  (5  ,  e~lt  I 
.  READ! 5  ,TH2*V0) ■ 

R E AD  (  5,  L  OS .» cs  ) 

CALCUtATIC.V  CE  PRELIMINARY  CONSTANTS 

Fl*3. 1415927 

G*32.2 

PR  A A=PR  H  20 

V  IS  CCA=VSCE2  0 
THC0NAsKN2  J 

PR  AN=PR  STM 
VI  <CCN  =  VSC',tm 
ThCCNN=k  S  T V 
X  ASUM  E=  J  .0  5 
EC  L=. 01 5 
GC  =  22.2 

V  AL  V  E0=  0  P°nT 
K  SCO  3  h  SO  rlO 
KSOF=KS  03C 
l«SDE*K3UF 
KSEF=KSOF 
CM- .11 
CPV=CI“ 

CP  Y=CM 
CPX*  c* 

CPZ=CM 

CPW=CM 

K3GH  =  KS  OCO 

KS JL  =  K  SO  f'S 

KSJK=NSJL 

PSH  L»KSJL 

ENTR  C0  =  E  STHD 

ENTFGH*  ENTOC 

BENDCDsAcNOHO 

AENOCh= OENJNA 

EX ITCD=l .0 

EXI  TGH  =  1.  0 

PS  AT  =  PCRUM  +  14.7 

HSAT  =  EXP  <0.26452 *ALCC( PSA'  )  *4  .46  70  3  J 

HFG=S22.  15-C.4G516*P  St  T>  l.  7  L  7E-04»  d<>  t'«*2. 0-4. 2i9E-03*°StT**3 .C 
FV  =  F  S  A  T  >  E  F  0 
H  M  M  0  *H  V 
FF  =  hSt  T 

TSAT«EXP(.22Ibl*ALCu(0SAT  )  >4.7  7123  ) 

MFCQt*Jl  IF  lU/3oGC.  J 
MACCO-  A  I-fFLT/  260 w .3 
TRBC*T„ASSC 
TVVO=TSCRN 

ThVO=T$PhTS 

TQRO-TkRC 
TSSO*T0ASSH 
TPPOaTS  AT 
TTTQ*TGA  S'>3 
TU'JO*T  0  A  S  EC 
TTLG»(TT10>TUU01 /2 .0 
TAMRa  00 . £ 

“MMO*T  OT *T  •> /3600  .0 

MMF6  .  000 1 

PHMI  n=  SRHST''/J600.0 

TNNO«S*CT 

TPPC=»0  SHOT 

PNN0»SFA3>14.7 

FPPO*FOPLH 

MMM  I  laJSHSTM/360C.0 

*NNO  **"'411 

TAAO*EC  CM  T 


C0N02 1  10 

C0N02120 

C5NC2130 

C3N0214J 

CONC^  150 

C0NC2  160 

C0N02170 

C3NC21 eO 

CON G2 l 90 

C0NC2200 

COM  0^210 

CON 02220 

CON02230 

CONC2240 

C0N022  50 

C0NC2260 

CONC227J 

C0N02230 

CON  C22SO 

C  CN02  300 

COMO  2  2  10 

C0NC222O 

CON 023  30 

CJNC234C 

CON  C2  3  50 

C0NC2J60 

CON  C 23  70 

C0N02360 

CON02390 

CON  0  2400 

C0N32410 

C0.NC2420 

CON  C2  43  0 

CONO  2440 

C0NC245C 

CCN0246U 

CDN0  24  70 

C0NC24c0 

C3N02490 

CON  C2  5  CO 

CDNC2510 

C0NC2520 

CON  C  2  5  30 

CCN0254J 

C0NJ2550 

CON  3  2  5cO 

C0N02570 

CO  N  C  2  5  e  0 

C0NC2550 

C0NC2600 

CON  C  26  10 

C0N02620 

CONC2620 

CON02640 

C0N02650 

CON  C266C 

C0NC2570 

CON 02  6  SO 

CON  C  26  50 

C0N02700 

CON02710 

CON  0  27  20 

C0NC2730 

CON C2 7 40 

C0MC2750 

C0NQ2760 

C0NC2770 

C0N02730 

CON  C 2  7 SO 

C0NC2aCQ 
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FILE:  CCFSTA7T  FCRTRAN  »l 


NAVAL  POSTGR AOUATE  SCHOOL 


C 

c 

c 

r. ' 

c 

c 


TRRO*EC  O^OT 

V8e0=.0160  048  8-.C000  0  20146*T3  AC*.  000 0000 36 5 U* T9  8 U** 2.0 
£-8.  142E-  il*TAo0**  J.O  ♦  !  .4GH  16-1 J*T  090**4  .0-1  . 14  3  E-  1 6  * 

CTBRO  **5  .0*o.u34E- 2G*T«9u**fc.C 
CTCTO*h3FVCL*LOOO.O/ooOO.O 
.  HNNO*HSHCUT 

CFUR *6=  F FAS *1000.0/ 3C0Q.0 

ARAC*RHA$SC 

HP  P  C*HL)  SO  L  T 

LCF1=L  AVSC 

LOF2=L3cl*NTtJBSC 

CDF* OT  J8  SC  *NTU3SC/12.  0 

COE*CCF 

OEF*ODt 

ACF=°1*CT'JBSC**2.C*NTU9SC/(4.C*144.0) 

VC  10F  =  A  CF* l UF1 
MASSV*‘-1i  SSiC 
MASSV^ASSS  H 
LJLl  *LA  vM 
LJL2=MTIR4^*LJL 1 
0  JL=OTtja*-t>*NrtJiiM  t/ L^  .0 
DJK*0JL 
CKL*0J7 

AJL*FI*CT'ie*e**2  .0*TUB^B/  (4  .0  *144.0  1- 

VJLJL*A JL*L JLl 

MASSY*  "•aSS1' 3 

*AS$>=MAS$tC 

NASSZ  =  **<3iOS 

ZBENC1* F  IIRM-hhOO 

LDF*2etaCl 

ZOF*  Za  ENO 1 

Z  8  EN  C2  *  F  a  CR  “-FRT  ROM 

L JL*Z9£  NC2 

ZJL-ZBENC2 

ZCC=ZCF 

ZGH-Zjl 

CROC  F0=0  .0 

DRCJL0*0. J 

LAB=LTUREr 

CAe* CTuefc/  u.o 

AR  EAEC  *fsTL8£C*LAB*PI  *OAB*  NFASS  E  '' 

ONP«CTJHOS/  12.0 

IC0*1AVHC*NTU 6nJ 

IGH*LA  *0  3*Nrj30D 

OCC*  CTU  9  F  0  *'l  rUBHC/  12.0 

OGH*t;TuBCL*‘jruaoc/i2  ,o 

AC0*NTg3H  J*  (  »  I*  (  070*00/24.  01  **2.0  ) 

AGh*NTU  9  C  C  •  I3  l*<  Cr  l“ 00/ 24.u ) **2.0) 

VOLDFM=t  Pi  *  JSTNOM**  2.0/  (  4.  0*14  4.0)  >  *LS™  O’* 

VOL  HJ*( o I*JMTR 0  ***2.  G/<  4.0*144)  )  *LWTRQM 
FiA  AO  *H6C  IN 
H860*HEC  CUT 

RHOF*toj .d-0. J 178 l*TSAT>  1. 1 32E- C5*TSA  T**2. 0-6.  786E-08* TSAT**3. 0 

VF*1  .0/  ArtCF 

VFG*  5  24.  C/PSA  T-0.  I 

VV*V  F*VFC 

RHOV*l . O/VV 

RHOFG*  RHOF-PH'J  V 

PHQM  *0*3  FCV 

RHONNO** SrlOLT 

RH0PP0*3 CSOOT 

HA AO 

M8B  C*M  A  A  C 

CALCULATE  THE  TCTAL  MASS  CF  FLUE  GAS  IN  THE  FURNACE 
MASSCR*RFLUE*FURV0L 

CALCULATE  THE  TOTAL  ENERGT  ENTERING  BOILER 


CCK02810 
C0NC2B20 
C0NC2E  30 
C0NO28  40 
CON  C2  6  50 
CON  02860 
C0NC2a70 
C3NC2880 
C0N02840 
CON02900 
C0N029  10 
C0NC2920 
CON  C  29  30 
C0NC2940 
C0NG29  50 
C0NC2960 
C0N029  70 
CON  C  29  60 
C0NC2990 
CON03000 
C0NC3C10 
CJN 030 20 
C0NC3O3O 
CONC3C40 
CON03050 
CON03060 
CON  030  70 
C3N03030 
C0NC3090 
C0N031CJ 
C0N031  10 
C0NC312C 
CON  03  I  30 
CON03140 
CONC3150 
CON  C  3 1 60 
CON03170 
CON  0  3  1 30 
C0N01190 
C  ONOJ  200 
CONC32  10 
C0NJJ22J 
C0NC323C 
CON  03240 
C0NCJ25J 
C0NC326J 
C0N03270 
C0N03280 
CON03290 
C0NO33OO 
C0NC33 10 
CON  C  3  3  20 
CC 7033  30 
C0NC334C 
CON  033 50 
C0NC33  60 
CON  033  70 
C0N03380 
CONC33  SO 
CON  03400 
C0NC34 10 
C0NC3420 
C0N03430 
C  070  3** 40 
C0NC3450 
C0N034&0 
C0NC34 70 
CON  034  60 
C0N03490 
C0NC35C0 
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FILE:  CONST ANT  FHqTRAN  0  I 


NAVAL  POSTGRADUATE  SCHCCL 


QQC*CTGTC*FURVCL 

CALC'lLATt  "L'llPEC"  LOW  £R  FEAT  IN  G  VALUE  OF  FUEL/AIR  Mix 
FHV*  CCO./MFOQO  • 

CALCULATE  ENERGY  ABSORBED  BY  RADIATION  IN  F'JRNACE-C 
C1C*CFLFA3*ARAC 
C2  0*C10 

CALCLLATE  THE  STEFAN-BCLTZFAN  CONSTANT  MULTIPLIEC  by  tfe 
SCREEN  AREA 

SIGMAAsQ  10/1  (  TOR  C*46C.  0)  *4  4.  0-(TVV0+46J.  0**4.  0) 

CALCULATE  THE  TOTAL  FLUE  GAS  MASS  FLOW  RATE  INTO  80ILE0 

PCCO  =  PFCCO**ACGO 
MRRU=MOGC 
MS  S0=M  R  RC 
FT  10  =  P  S  SC 
MUUC=HT1C 

COMPUTE  THE  ENERGY  TRANSFERRED  TO  the  WATER  (LOWER)  DRUM 
BY  THE  DESUPERHEATER 

O9O=NNNO*(HNNO-HPP0) 

Q990=C9J 

CALCULATE  THE  ENERGY  TRANSFERRED  TO  tH£  MAIN  BANK  RISERS 

C6C=FP*0* (HMMQ-hEEO ) -C20-0990 
Q  5  C*  Q  6  C 

COMPUTE  THE  PA7N  BANK  HEAT  TRANSFER  CO  EFF  I  C I  ENT-TJ  BE  MET  A  L 
TO  STEAM  SIUE3 

KY*AREAM3/1.  782EC6 

CCPFJTE  THE  main  BANK  TUBE  METAL  TEMPERATURE 

TYY0=(060/(KY*oSAT**(4.C/3  .0)  )  )**(  i  .  0 /3.  0)  ♦- T S AT 

COMPUTE  THE  MAIN  BANK  HEAT  TRANSFER  COEFF I C I EN T- FLUE  GAS 
TO  TUBE  metal 

TSTG=<  T $  SO  +  TTTO ) U.  G 

KST*C50/  ('■>SS0**0  .6*  ( T  STU-TYYG  ) ) 

COMPUTE  THE  SCREEN  BANK  HEAT  TRANSFER  COEFFICIENT-TUBE 
PETAL  TO  STEAM  SICES 

KV*C10/(PSAT**(4.0/3  .0)*<  TVVO-TSA  T)  **3.  0) 

COMPUTE  THE  HEAT  TRANSFER  TO  THE  SUPERHEATER 

04C*«MMO*(HNNO-HPPJ ) 

030*040 

COMPUTE  THE  SPECIFIC  HEAT  OF  THE  STEAM 
IN  TFE  SUP £3  PEA  TER 

CMN*040/(MMMO*l  TNNC-TMMO) ) 

COMPUTE  THE  SUPERHEATER  HEAT  TRANSFER  COE FF I C  I  ENT- FU J 5 
GAS  TO  TUBE  METAL-C 

TRSC*<TRROTSSO>  /2.C 


CON  035  10 
CONO  35  20 
CON  C 3  5  3 C 
C0NO354O 
COfJCi  5  50 
CON  C  3  540 
CON  03  5  70 
CONO 35  P 0 
CON C3 5  SC 
C0N03o00 
C0NC3610 
CON  C  34  20 
CGNu3630 
C0N0364O 
CONO 30 50 
C0N03660 
CONC  34  70 
C  ON C  3 4  SO 
CONO 3690 
CON  C37CC 
CON  03  7  10 
CONGA  720 
CON  C  3  7  30 
C0NQ374O 
CONO 3 7 50 
C  ON  03  7tO 
CON03770 
CON C3 7  SO 
CON  C3  7  JO 
CONO  33  00 
CON  C  3 1 1 0 
CONO 3320 
CONC3BJO 
CDMC3E40 
CQN033  50 
C0N03B60 
C0N036 70 
CON038M0 
CON C 38  SO 
C3NC39C0 
CONC3910 
CONC3S2C 
C0N03930 

conc39*o 

CON  C  39  50 
CC NO 3 9 60 
CON03S  70 
CON  C  39  BO 
CON03990 
C0NC4C00 
C0NC4O1O 
CON04023 
CON  C4C30 
C0N04040 
CON  04  05  0 
C0M040o0 
CON04070 
CONC4CcO 
CON  040  90 
C0N04100 
CON  C4  110 
CDNC4120 
C0N041 30 
C0NC4  14C 
CDN04  150 
C0MC414C 
CON  04  1  70 
C0NQ4180 
C0NC41SC 
CON  04200 
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nnoo  Onn  non  nno  noon  onnn  nnn  non  r»°oo  nnn  non  nno  nfin  nn^n 


NAVAL  RQSTGRACUATE  SCHOOL 


:  CONSTANT  fortran  PI 

KR S*  030/  <MBR0**0.  6*<  TRSa-TxwO)  J 

COMPUTE  the  heat  transfer  to  THE  ECONOMIZER 

cac*nto* (heao-HAAT j 

0  7C=080 

COMPLTE  THE  SPECIFIC  HEAr  CF  THE  FEEDWATER 
IN  THE  ECONOMIZER 

CA8*08C/I *4AC*( Ta80-TAA0) ) 

CCMFLTE  THE  SPECIFIC  HEAT  OF  FLUE  GAS  IN  THE  FURNACE 
COR*  (COO -010  J  /  (mrrj^iTORO-  T4MB  )) 

COMPUTE  THE  SPECIFIC  HEAT  OF  FLLE  GAS  IN  THE  SUPERHEATER 
CRS=030/(MRRC* ( TRRC-TSSOI » 

COMPUTE  THE  SPECIFIC  HEAT  OF  FLUE  GAS  IN  THE  MAIN  BANK 
CST  =  050/ ( MS  SO  *(  TSSJ-TTTQ  )) 

COMPUTE  THE  SPECIFIC  HcAT  OF  FLUE  GAS  IN  THE  ECONOMIZER 
CTL»070/t »TT0*(TTT3-tUU0I ) 

CCHFUTE  THE  ECCNC-'IZEP  HEAT  TRANSFER  COEFFICIENT-TUBE 
METAL  TO  LI CUIu 

KX=(  U023*TftC0NA|/CAE)*(4  .0/  I  P  1*0 A3*V  IS  COA*NTUBEC)  J**C.8 
e#PRAA**C.A*A:’  EAFC 

CCMFLTE  Trie  LMT  C  FOR  T HE  ECONOMIZER 

LMTOAe*08I/ (KX*M* 40**0.81 

COMPUTE  THE  ECONOMIZER  TUBE  METAL  TEMPERATURE 

E  XP0EC*E  J>°  (  (  T33J-TAAC)  /LMT0A8) 

TXX0*IT AAO-TeB^*EX30  EC )/<  I  .o-EXPOCC  ) 

COMPLTE  THE  ECONOMIZER  HEAT  TRANSFER  COtH FI C I EN T-FLUE 
CAS  TO  -U8E  METAL 

KTU*Q70/(  M  TT 0  *(  TTUG-TXXO)  ) 

COMPUTE  THE  DE SUPERHEAT  £R  HEAT  TRANSFER  CCEFF  I  C  I  ENT-T'J  B E 
METAL  TO  STEAM 

KZ*  I  (  .02JMTHC0N.N  ) /JNP) *( 4.  0/1  Pl*ONP*  VI  SCQN*NTU3DS)  J 
C**0.8*°R  AN**0 . 3* ARE A CS 

COMPUTE  THE  DESUPERHEATER  LMTO 

LMTONP*QSO/(KZ*MNNC**C.3) 

CCMFUTE  THE  SPECIFIC  HEAT  OF  STEAM  IN  THE  D£SU° ERHEATER 

CNP»C90/ (MNNO*ITNNJ-TPPO )) 

COMPUTE  OESLPERHEATER  TUBE  METAL  TEMPERATURE 

E  XPQO  S  =  E *P{ I TNNG-TPPCI  7L  MT  CNP  J 
T  ZZO*  (TNNJ-  tXP03S*T°°w)/(  I  .U-EXPODS) 


) OS  *T°  n w I / ( I .0-E  XPOO  S) 


COMPUTE  THE  05SLPERHE ATER'  HEAT  TRANSFER  CCEFF I C  I  ENT  - 
WATER  ORUM  LI'JUIO  TO  TUBE  METAL 


CON04210 
CON C 42 20 
C0NQ4230 
C0NGh240 
C 3 N C42  50 
C0N0426J 
C JN042  70 
COM  042  80 
C0NG4290 
C3NC43C0 
C0NC4310 
C0N04023 
CONC433C 
CQN0434J 
C0NC4  350 
C0NC436Q 
CCN0*»37  J 
C0NC43E0 
CON  C4390 
CQN044-00 
COM  044  10 
C3N04420 
C0NC443C 
C  3 NO 4 4 40 
CCN04450 
CONC4460 
C0NC4470 
C0NC448J 
C0NC449C 
CON  045  00 
CONC4510 
C0NC4520 
C0NJ4530 
C0NG4540 
CON  045  50 
C0N0456J 
C0NC4570 
C0NC45F0 
CONC^bgo 
C0NC46CC 
COM  046  10 
CJNC4620 
CON C 46 20 
CONG4640 
C0NC465C 
CON04660 
CON04670 
CON  C46  60 
CONC46SO 
CON047  00 
CON047  10 
C0N04720 
CONO  4  7  30 
CON  C4  740 
C0N0475J 
C0NC476C 
C0M04770 
C0NC4780 
C0NC47S0 
CON04800 
C0NC48U 
CON  043  £0 
C0N0483C 
C0NC4840 
C0NC485C 
C0NC4860 
CON  C4E7C 
CQN04860 
C0NC4890 
C0NC49CO 
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FILE 


CONSTANT  f  CAT  RAN  PI  NAVAL  POSTGRADUATE  SCHOOL 

KHJ=C99C/(  T2ZC-TSAT) 

COMPETE  The  S  J  P£  PhE  AT  ER  L»'T0 

LMTCNN=(TNN0-TMM0  )/  (  ALOG  {  t  TWWG-T**mO  )  / (  TM V.O-TNN  C >  >  ) 

CALCULATE  THE  SUPERHEATER  HEAT  TRANSFER  CCEFFI Cl ENT-TU  E 
METAL  TO  ST  t  AM 

KW*Q40/(MMM0**0. S*LmTQNN) 

CALCULATE  THE  THPJTTLE  VALVE  FLCw  CCEFFIC IENT 

K0N4-MMMII i /( VALVE)*PNNO  ) 

CALCULATE  THE  SUPERHEATER  nUTLET  DENSITY 

K0M3=(  PNvtO^RNMJ)  /(  IhCMMO  +  RHONNC) 

PHCNNO*  (  0  n  C  M  “•' )  ♦  »KNNO  )/2  .0 

K0N1«( ( PMNL-FNN0)’»AnCMN0) /MMMO**2.0 

FOP  THE  DOWNCOMERS 


CCKPUT&  ThE  FRICTION  PACTCRS 
ANO  THE  RISERS 

FC  D-*  1 . 0  / (  1  •  7 4-2.  0*4LC3L0  KSCC)  ) 
FGH»  1.0/  (  1.  74-2.  C*\L0G10<  <  SGH  )  ) 
F0E»1.0/ <1.74-2. O^ALEGIO  <<  3  U  6  )  ) 
FEF=  l.  0/1  L.  74-2.  Ca  4L  OG 1 0<  K. SfcF  )  ) 

F  JK*  1 .0/  (1 . 7+-2.U  *  iL  TG1JI  nSJM  ) 
FKL*l. 0/(1. 74-2. G*4lCG10 (K.SKL) ) 


START  ITERATION  TO  BALANCE  CIRCULATION  LOOPS 


COMPUTE  INITIAL  VALUE  of  RISER  OUTLET  MASS  FLCW  RATE 

MFFC=0  10/1  XA  SUME  *HFG  ) 

MLL0=C5O/(X  4SUME*PFG) 

CALCULATE  the  initial  downcomer  ENTHALPY 

71  HCOC*(  (  MFF  0+MLL0-M3dJ>*HF+NSe0*h88040  .0  *  (MFFO  MMLT.0  ) 

6  *HV  )/(  MFFU+iMLLO) 

HGHO  *hCC0 

COMPUTE  THE  INITIAL  ORUM  ENTHALPY 
FORUMG«HCOO 

COMPUTE  MAIN  3  A  NX  RISER  I NLET  ENTHALPY 

MhFO  *  ML  LO 

H J JC*HGH  C>  090 /MHHO 

CALCULATE  THE  RISER  INLET  DENSITY 

RHCDOO*RHDF 
RH  C  J  J0  =  PHCF 

COMPUTE  THE  COnNCOMER  OENSITY 

VCCC* ( (  MFF  C+MLL 0— M3B C)  *VF>  MB  30*V3 BO) 

6/  (MFPO^'U.lO  ) 

RHOCOG-=l  .9/VCOO 
RH0GF0®5  HOCOO 

£  \  CALCULATE  THE  RISER  OUTLET  CUALITY 

C 


CO  N  0  4  9  1 J 
CON0492J 
CON  C4  4  20 
CON04940 
CONC4S5C 
CON  04960 
CCNC4R70 
C0NC49HO 
C0NC4990 
CONOSOUO 
C3NCSCIC 
C0MC5G20 
C0NC5G30 
CON  C  SC 40 
CON05O5J 
C0N05060 
CON  0  50  70 
C0N0508J 
CON  C5C90 
CON C 5 1  CO 
CO  NO  SI  10 
C0NC512C 
CD  N05130 
CONC5140 
CON  C  5150 
C0N0516G 
CONC517C 
CJN05L30 
CON05I90 
CONC52CC 
CON 052 10 
CONC522C 
C0NC522G 
CON05240 
C0NC5250 
CON  0  5260 
CON05270 
C3NC5260 
C0N05290 
C0NC53  00 
CDNC53  10 
C0M05320 
C  ONO  5330 
C0NC5240 
C0NJ5350 
CON  C  5  2  60 
C0N0537O 
CON05330 
CON05290 
CON05400 
C0flC54  10 
CON  0  54  20 
C0NO543O 
CONC5440 
CDN05450 
C0N05460 
CQNC54  70 
C  ON  054  60 
CQN05490 
CONC55CO 
CON  055 10 
CON  C5520 
C0NC5520 
C0N05540 
CON0555O 
C0N05560 
C0NO55  70 
CONC55EO 
CON C5 5  SO 
C0NQ560O 


UO'-J  WOO  UoW  Woo  WOO  WOO  OwW  OWOW  wowvj 


FILE :  CONSTANT  FORTRAN  PI 


NAVAL  POST GRADU ATE  SCHOOL 


HO  CO  =HC CO 

XFF0*(  JIOMFFO*  t  hF-H COO  )  )  /  l‘xFF0*HFG) 

XLL0= (  05U*'1LL  J*(  HF-HJJO)  )  /  <  ILL  C  -MFC  ) 

CALCULATE  THE  RISER  OUTLET  JENS  ITY 

RHCFFC=SFCF-XFFO*RI-Hpa 
RHOLL  C=PHOF-XLL  0=*R‘H0  FO 

CALCULATE  THE  NCN30I LING  LENGTH  PF  THE  RISERS 

2BENC1=  FNOEM-HHCP 
LDF  =  Z6EHCl 

LOEO=LOF  «AM  AX  K ( HF-HCDCJ ,  0.0)/{ (HOD  0+ XFF C*rtFG) -HD00) 

2BE  FC2SH N C  HWT  c  CM 

LJL*ZP£NC2 

LJKO*LJL*AMAXI(  (HF-HJJO)  ,0.0)  /(  (  HJ  J  0-t-XLL  0*HF  G  )  -H  J  J  0) 

CALCULATE  SOILING  LENGTH  CF  RISERS 

lEFO-LQF-LUEO 
LKLQ=L JL-LJXC 

CALCULATE  THE  AVERAGE  DENSITY  INThE  SCREEN  RISER 

PFOCFO*  (1  .0/LCF  )  *  (LE  FO/  (  (X  FFO  )  ♦VFG)*ALr»G(  (  I  XFF  C 
£ )  /  VF  )  *  VF  3* 1 . 0  )«-  R  HC  020*  LOE  0  ) 


CALCULATE  THE  AVERAGE  DENSITY  IN  THE  MAIN  94MX  RISER 

RHCJ  L0= (  1.0/L JL J  *<LKLC/(  < XLLO ) * VFG ) *AL0G ( ( XLLO 
£* Vf G)  /  VF  +  1 «  3  1  ♦KriC » „0  *  L  JKO  ) 

CALCULATE  THE  EFFECTIVE  HEIGHT  OF  THE  RISERS 

ZDE0=Z0F-LEF0 

2JXO=ZJL-LKLO 

ZEFO=ZCF-ZDEC 

ZKLO=ZJL-ZJXO 

COMPUTE  THE  TVO  PHASE  FLO*  MJ  LT  I  FL I  C  ATI  ON  FACTORS 

PG  PA VE*29.  794*  X  F  F  3**2  .0-6.5C66  *XFFO  +  .9776 
RGR  AVK*  2  A.  794*XLLC**2.  0-6.  5o66*XLL0+.  9776 
PACLE*15.456i-*XFFJ**2.0*-1o.49‘,4*XFF0—  .00007 
RACL  J*l  5.4564*XLL0**2.O*l3.4944*XLL0-.000O7 
RFR  ICE*- J*.0S22*XF=0**2.0*23.  71o4*  XFF  0+ .  3  734 
RFR  l  CK*-34. 0  322*  XLLO  **2.0 +  23.7  lfc4*XLL3+.373x, 

CALCULATE  SFC3N0  APPROXIMATION  CF  MASS  FLO*  RATE  AT 
EXIT  OH  SCREE  I  S  IS  ERS  M 


PHDEEO’RHROOO 

PFF  CC* ( ( PHCCGO*G*Z:C-G*ZCEO*( ( 5 HOD DO ♦ FHDF  EO  J /2 . 0 ) 

£— G*Z  EFO*R  FOES  0*R  GRAVE  )  /(  (  r  C  0  *  LCD/JCJ  +  E‘JTCCJ 
£*eE'.  CCC*E<  lrCC  )/  (  c.O  <ACD**^  ,0*'’HO:  3  0)*(  (=H  'EEO-R  HODDC  ) 
£  /(RHPEEO«RHC 300*A0F“*2.0 J ) + < 4 . 0* FOE* L DEO  *2  .0 
-  6)/(2.0*DOF*C!HTEEC«-r  HDOO  0)  *ADF  **2.  C)  *'=  AC  LE  / 

£  (RHCEtO*ACF**2 . J  )  ♦ (4  .0 *F E F *L  E F3 » F=  ICE )/ 

£( 2. 0*D0F*Rh0cE0*mDF**2.0) II **3.3 

CALCULATE  SEC3NC  AP  FR  CX  IM1TI3N  CF  M  AS  S  FLOW  RATE 
AT  EXIT  OF  MAIN  3ANX  RISERS 

RHOKKO*PHOJJO 

MLLOO*  (  (  R  H’,)GH,0*G*Z  3H-3  *Z  JX  0*(  (RH0JJ3+°HDKK0)  /2.0) 

£-G*  ZK  LD  *  AhUKKJ *  R  «R  AV  ’<  )/  (  (  F  CH*L  CF/ I3GH  +  E'I  7  =  GH 
s  E*3EN0GH*E  XI TGH) /( 2. 0*4 Grt* *  2 . 3* RH CGH 0 ) + (  ( Fri-JXXD -RHCJ JO  ) 
.  £/  (RHCKXJ  *RH9  JJ 0*4 JL  **2.0  )  I  ♦(  4 . 0 *F  JK *L JX  C «  2.  0 
£ )  /  C2. C*0 JL* IRnJMO  *R  HO J JO ) *AJL**2.0) ♦RACLJ/ 


C0NC5EI0 

CON056  20 

C0ND56  30 

CON  C564C 

C0NC5650 

C3N05660 

CONCEfc  70 

C0N056  80 

C0NG56SJ 

CON  CS  7  OJ 

C0N057I 0 

C0NC5720 

C0N05720 

C0N05740 

C3NC575C 

CDN057UO 

C0NC5770 

CDNC5760 

C0N05790 

CON05300 

CON  058  10 

C0N05  620 

C3NC58  20 

C0NC6F4G 

C3N053  50 

C0NC566C 

COM  058 70 

C0KC583O 

CON  C5ES0 

CDN05900 

C0NC59I0 

CONG 59  20 

CCN05930 

C0NC5940 

C0NC5S50 

C0N0396 0 

C0NC5970 

CON  0  59*0 

C0NC5990 

CON  C6GC0 

C0N060 10 

CON06C20 

C0ND6030 

CGN06040 

CQN06C50 

C0NC6C60 

CCN060  70 

COMC6C6C 

CON  06090 

C0NC61 00 

CO  N  C  6 1 1 0 

C0NC6120 

C0N061 30 

*0NCfcl4C 

C0N06150 

C0NC6  160 

C0N06170 

CON06 183 

CONCtlSC 

CON  06  2  00 

C0NC62 10 

CON Ct 2 20 

C0N062  30 

CON06240 

COM  06  2  50 

C0N06260 

COMC627C 

CON 06 2 80 

C3NC8290 

CON  C62C0 
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l 

c 


6CQ 

601 


£IRH0KKC*4JL*"2.G)  «■<  4 .  C*F  KL»  L  K  L0*  RF  R I  CK)  / 

C0N06310 

£12.0 *CJK *RhlK.1sJ»AJC  *»2.0 ) J  )**0.5 

C0N0632 J 
CON 06320 

COMPARE  PREVIOUS  APPROXIMATION  FCR  RISER 

MASS 

C0N06340 

FLO*  cAT£  TC  CURRENT,  IF  WITHIN  E=R0R  CRITERIA  CONTINUE, 

C0NC625J 

IF  not,  update  ANO  t HEN  REITERATE- 

C  ON  06  360 

• 

C0N06370 

CHECK*0.0 

CJNC63EG 

IF (ARSl  IFFu-MFFCC)  .LI. .01 I  GO  TC  52 

CON  06  390 

MFFQ*  <MFFOJ--'FFO» /2.0+MFFO 

CONC64  00 

CHECK*! .0 

C0NC6410 

1  IF(  A8SI  MLLC-RLLOC) .LT..01 ) GO  TC  54 

CO  NOo  4  20 

MLLO*  IMLL  J0-MLL0  >U  .0*‘U.LJ 

COM 06430 

GO  TO  71 

C0N0o440 

•  IFICFECX.EO .l.OIGG  TO  71 

C0N06450 

CONC6460 

COMPLTE  INITIAL  MASS  OF  LIQUID.  IN  DRUM 

CON06470 

C0NC6480 

DMASlO*  (  VCLQRMM.RFCCCO  J/2  .0 

C0NC64SC 
CON  065  JO 

COMPOTE  ThE  11ITIAL  ORUM  "LIQUID"  VOLUME 

C0N06510 

C0N06520 

OMOV  C=  VOLORM /2 . 0 

CGN06530 

m 

C0NU6540 

CCMPLTE  ThE  INITIAL  ENERGY  STORED  IN  DRUM 

L  IQUU 

CON  Ce 5  50 
C0N0&560 

CM Ch  10*  C K  AS  LO  *H C RU MO 

C0NC6570 

C0NC65S0 

EQUATE  INITIAL  FLOW  RATES 

C0N06590 

C0NC66CC 

MCOO*MFFC 

CO  N  0 1 6  1 0 

MGhO=MLLO 

C0N06620 

COMPLTE  INITIAL  MASS  CF  STEAM  IN  STEAM 

CONC6620 

CON06640 

CRUM 

C0NC6650 

C0.N36O60 

OSTMO*VOLORM*RHOV/2.  c 

CON06670 
C0N066  60 

COMPLTE  HXFFO  A  NO  HXLLO 

CON  066  SO 
C0N06700 

HXFF0*FF*FFC»XFF0/2.0 

C0NC671C 

HXLLC*HF«-HFG*XULC/2.  C 

CON 06  7 20 

WRITE1S  ,600  ) 

C0NCo730 

CONC674Q 

FORMAT!  1H  l) 

C0NC6750 

WR IT  E(o,  IN CON  I  J 

C0N067 e>3 

WRIT  E (6 ,601 ) 

CON  Cfc 7  70 

FORM  >T(  I.-I  O) 

CON067bO 

WRITEI6,  I  ICON21 

C0NC6790 

WRI  TE! 6  ,601  I 

CDNC68CQ 

WR  ITE ( o,  I  ICON  3) 

C0N068 10 

WRIT  E (6  ,60 1  ) 

CON  C6E2C 

WRI  TE  (  6  ,C0  jST1) 

C0MC683C 

WRITE(o,6UI) 

C0NC6840 

WRITEI6.CCNST2) 

CJNC6650 

WR ITEI 6, 6CI J 

CO  N 06  860 

WRITE(6,C0nS”3) 

C0N068  70 

V>RI  TE  I  6  ,601 1 

CON  06880 

WR  IT  6(  o ,  C -IN  ST4 1 

C0N0689D 

WRITE  (a  ,60  I  J 

CON C6 SCO 

WRITE!  O.CCN3T51 

CON  06  S 10 

WRITE!  o,60  1) 

CON06920 

WRI TE (6 , 600 ) 

CON  C6S3C 

WR ITE! 6  »OuTQ  UT 1 

CO NO 69 40 

WRITE17,  IN  CON  I ) 

CON06950 

WRITE!  7,1  NCCN2) 

C0NC6S6 J 

WR  ITEI  7,  INC  AN  3) 

CO  NO  6  9  70 

WR IT E  1  7  ,  C  Ci',S T  l  i 

CONCfcSEC 

WRITE! 7.CCNST2I 

C0MC6990 

W  R  ITE  !7,  CONST  j ) 

CQN07000 
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APPENDIX  B 


C6/3C /7<-  21.26.19 


FI  IE  :  CS  FF  FC9*  RAN  ?1 


NAVAL  RO STGR 4JLATE  SCHOOL 


//HOPEFULL  JD5  (  2473,C567,NE  61  »»9)  ,  '  WALK.E  R  c  ‘*C  1319*  ,TI«E  =  2 
//  EXEC  CS^°XV 
//>.  S  V  S 1 N  DO  * 


*  THIS  0  T  Y°  E  R 01 LE  R  MODEL  IS  WRITTEN  IN  CSFP-III 

*  CC  **  FJ  )T  E  R  LANCJACc.-  IT  is  4  'n--'LINE«P  M  ID  EL 

*  WHICH  LT1LIZEE  A  SET  OF  initial  ccnditiqns  anc 

*  CT  HEC'  DATA  PREPARED  9  V  AN  I N I  T  U  l  CONDITION  FORTRAN 

*  FCCEL  WF1C8  PRECEDE*  THIS  ONE.  T ale  w  )DcL  IS  GOOD 

*  FOR  SM  ALL  TFAhSIENTS  AF.  DU  NO  ThE  CPtRATINO  POINT 

*  SPECIFIED  3 Y  the  INITIAL  CDNDITIJN  PROGRAM. 

INITIAL 

THNST  DRA*1=C.  C»DP.  A42  =  6.0*0FA4  3  =  9.3,  £33*4=0.)  , CRA45=0 .0 , CR  A4b=0 .0*  ... 

D  R  74  7=0.0*  C  -  A4d  =  *  .J  *  OR  -  =  1=  C.3,D-''~*32=6.*,DRA;3=0.  0,DRAa4=C.C,... 

OR  A3  5  =  l.  C  ,L'F  A3  6=  d.  0,  OF  73  7  =0 .0  ,  CR  433=0  .0 
CONST  0't£=  l.  C,  ... 

TRP,0=  2563  .0000  ,TWO=  cOO.OOOCO  ,Tv»w3=  1033.  OOCC  ,TYYC  =  ... 

515.42745  ,T220=  600.99976  ,TXXO=  360.36450 

CCNST  TWO* 2.0*  ... 

0  F  J  H  0  =  2  154667. 0  ,C<ASr.D  =  5415.3164  , haAJ=  45.00C0J0  ,TA40=... 

2  54.  C  C  C  C  c  ,**FOOJ  =  3.307*999  ,MACC9  =  55.  398330  ,  MNNO*  .  .  . 

7.7500000  ,  *M*J  =  *3.000000  ,VALVEO=  0.50949939  ,3000  =  ... 

55.2C6375  ,CMN  =  O.o5523023  ,CA3  =  1  .0/6i9J0,CoV=0.1J999995»  ... 

0^=  0.10999995  *fPx=  C.  10999995  ,C  PZ=  0.10999995  ,CPw  = 

0.  10999905, GFJnLJ=39ol62o.J,  C'*  AS LO =3002  .20 JO, DM JVO=  163.65698 
CONST  Th°  EE  =2. C,  ... 

DRO  CFO=  0.0  *  )c  °  JLO=  0.0  ,MM-1E=  0.  449  9  9979E-02  *  .  .  . 

MGhO=  o62. 70.925  ,M9°0=  45.JJ0J00  ,LST'1D''  = 

U. 669996  *  9 S 1MDM=  6C.C0090Q  , MC 00=560. 02 734, CLCV 0=163 .66 

CCNST  F  ”'1 1  R  =  4  .0,  ... 

N  Sj  l=  C.5S5SS591E-D4,l*SJF  =  0 .29999996 E-04 , <S GH=  0  .4  5999  999  E-05  ,  K.S  C  D=  .  .  . 
G.49599995E-05.KS JC=  0. 2 9 S59996C -04 , NSE F=  0. 299999 96 t-04 ,  ... 

KSjl\  =  0.3S595  R91E-0*,  ... 

NSNL=  C. £99599915-0*. AJu=  10.678426  ,DJL=  193.57166  ,LJL2=... 

26420.000  , L  JL 1=  lC.CCJUoJ  ,DOF=  3 5.  89 747c , lDF= 1 5.  75 0 JOC  ,  . . . 

LDF 1=  12.7CCC0)  ,  A  CF=  *.06  46043  , L GH=  25.00O0J0  , A  G H= ... 

1.44  =  13  32,  SAMRl= 250,3  I STIM*  10.0 
CCNST  F  I V  E=5  .0,  ... 

DGH=  1.9183321  ,0C0=  4.6699991  ,LCC=  9*. 500000  ,LDEJ=  ... 

2.6166=15  ,LJNO=  C.  fcC  ECO  102  ,ACD=  2 .  95  47  77  3  ,  ‘-'A  S  S  V=  7951  .  3  9  34  .  .  . 

,MASSW=  10752.000  , v  AS  S  Y  =  2  =  192  .000  ,**AS.,X  =  2444 3 . 0 JO, W A SS2 =. . . 

lcl/.LCCC  ,Z3E9Cl  =  15. 75  OCOO  ,FCD=  0  .8  1 J  23  69  3E-0  1  ,  ... 

F Cc=  0  .9271*3465-01,  FEF=  C  .92714  =  465-  Ol,FGH  =  C.  8102365=£-C1 
CCNST  SIX =6.0,  ... 

Z9EWC2*  U.41CCCI  »8A  SS  JR  =  71.  395  ,VDLJL=  l  0  6 . 7342  6  ,  VC  L  C  .  .  . 

F=  51.620*68  ,  V  “L  JR  *'  =  227.JU96  ,V(DLHJ=  98.  1  74744 

CONST  SE  \EN  =  7 . 0  ,  ... 

CQ3=  J.jl5‘59al4  ,crS=  C.  3045367C  •  ,C  ST=  C.  29042  3J4  ,CNJ  = 

C. 26  02  6  064  ,CNP=  0  .533471  39  ,<NP=  13.506820  , K W= 2 . 095 7 7? 4 , . . . 

S  IG MA A=  c.  21  l22E12E-C5,lsR  S=  0.87139436  ,<rT=  ?. 2037563 

CCNST  E  IG  FT=  3  .0 ,  ... 

KTU  =  C.441266CC  , £N TP GF =  0 .99995 964 E-0 l , F NTPCO=  0.  9959 9 S64E- 0  1 ,  .  .  . 

3ENDC  ii-  C .  S599  4  Ft'.fc -  0 1 , E  X  I TGrt  =  1 . 000  JO  CO ,  C  <  1 TC  0=  .  .  . 

I  .0000000  ,  HA  A0=  22  1  .00000  .ND'1  1=  0.  1064  =  6  7oF-  C  l  ,..  . 

K.JN4=  .  1C777515  ,FJx.=  0.9819'  145E-01  ,F<L=  0  .9  3  l  It,  1*9  E-0 1 ,  .  .  . 

GC=  J2 . 1999  9 1  ,TAM8=  dC. CCCOCu  ,ZCD  =  15.  75  0000  ,2DF=  ... 

15.750000  , 2 JL=  12.410001  ,IJL=  12.410001  ,... 

3  END  GH  =  C.  1451SS9d,r,CN3  =  6*7.  714o0  ,«=32 .1 999 97  ,Q ST ''0=2*0  .59802 
CCNST  N1'E=9.0,  ... 

KX  =  1 0.  C  7  73  7  3  ,K2=  14.930669  , KV=  J . J2 7696 1 3E-C 5, . . . 

F  HV  *  1S3CC.  456  ,.ort=  35. 89  747c  ,DKL=  Ufa. 5716b  ,DEf=... 

35.697*76  ,CJK=  Ho  .5  7 1 66  ,r.  V=  0.  36683  50  2E-  02,  .  .  . 

Oft  A1C*C.C»0RA2  0  =  0. 0»XfF)  0.139503  00  E -01, ... 


CS'*3D0  10 
C  SWC0020 
CSMCCC20 
CSM00040 
C  S  **  C  CC5C 
CS-00060 
CS  "00070 
C  SWCCC60 
CSMuCOSO 
CS«00100 
C  SM  C  Cl  10 
CS"0J120 
CSMG0130 
CSM00140 
CS^'JOISO 
CSMCC16C 
CS" 00  1 70 
CSN00180 
CSMCC150 
CS  m  CC2  CO 
C  S  MOO 2 10 
C  SMCC220 
CS  **00230 
CSM00240 
CSMC0250 
CS  MC02  60 
,C  S**CC27C 
CS  M00280 
C  S  m CC2  90 
CS " CCD  CO 
,CS MOO  3 10 
CSMCC32C 
CSM00330 
CSMC0340 
C  SM  C  C  3  50 
C  >*-'00360 
CSM00370 
,ESM0O33O 
CS  mo 0390 
CSMCC400 
CS**J0410 
CSM00420 
C  S'-C0430 
CS  **00440 
CSMC0450 
CS '*00460 
CS  MJ04  70 
C  S'*CC460 
CS" 00490 
C  5MCC5CO 
CS**C05  10 
.  C  S  MO  05  20 
C  S**C  C530 
CS*‘00  540 
CS  M  C05  50 
C  S'*  CC560 
CS  HJ 05  70 
CSM00580 
CSMCC5SC 
CS  MQJ600 
C  SMC06  10 
CS " 006  20 
CS  “J 06 30 
CSMCC64C 
CS**C0c50 
C  S**CC660 
C  5**  CC6  70 
CSM0J680 
C  S  MC0690 
CSM00700 
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fjl£:  CSPP  FORTRAN  PI  NAVAL  POSTGRADUATE  SCHOOL 

*L LO=C.4tCCCCCCE-Cl#MFFC=S EC. 02734, MLL0=662. 70825 
CONST  TEN=10.0,... 

OR  A1  1=C. C, DR  A  12=0.1  ,DR  A  1  3  =  C.  C,DR£1  4=C.O,  D  R  115=0.0  ,  CR  Alo=o  .0 ,  ... 

04  A 1 7=0 • 0  ,0=  A2l=  0.0,3-  4  2  2=  C.  0,0  R  A  23  =  0.0,0  RA2  4  =  0.0,... 

FC L=. 015  ,CC0  =  63836. 345  ,  2GF  = 12  .410001  «... 

T>A2S=C.  C,0=  A 2 6=  0.  0,  JR;-2  7=0.0  ,riXLLC=  505.  81  ,rtXFF0=502. 34 
DYNAMIC 

*  INPUT  EQUATIONS 

VALVE*.  015*RA*'PU  0. 0)-. 015  *RAMP  (20.0  )*.5l 
«>NNs  MNNC 
*1A=MAA0 
TAA=  TA AO 

*  CCMFIJTE  Th£  TOTAL  FLUE  GAS  FLOW  RATE  INTO  BOILER 
KC*“CCO 

*  COMPUTE  THE  ENERGY  ENTERING  THE  eOI  LER 
CO  *QQO 

*  CONFUTE  the  ENERGY  TRANSFERRED  TO  THE  SCREEN  RISERS 

*  VIA  RADIATION 

Cl  =SI  G  FA  A*  IlTRR+460  .01**4.  O-ITVV  +460  .0  )  **4  .0  ) 

*  COMPUTE  THE  RATE  EQUATION  FOR  FURNACE  FLUE  GAS 

*  TEMPERATURE 

DT  RR=  (C0-C1-MRR*C0R*(TRR-T  AMg  J  )  /  ( *n  S  SQR*C  OR  I 

*  COMPUTE  THE  flue  gas  temperature 

TRR*INTGRLMRRC,DTRR) 

*  CC“Pre  TF£  •rE,-FE=ATJcE  OF  the  FUJE  gas 

*  LEAVING  THE  SJPtRHEATfcR 

PHIl=?.C**FH**C.4«CHS/KP5 

TSS=  I  TR°*(  PH  I  I- 1.01 +2.  C*TV>  «) /I  PH  1 141. 01 

*  CC*PUTE  THE  SU PE Rh.E AT Eft  ENERGY  TRANSFER 

*  FLUE  GAS  TO  T'JR £  METAL 

Q3*MRR  *C  R  5*  (  1RR-TSS) 

*  CC v  FU" E  Thfc  '31'  BANK  ENERGY  TR ANS FER 

*  flue  gas  tj  tub c  metal 

Q5»“S:*C  it*  nSJ-TTT) 

*  compute  the  temperature  of  the  flue  gas 

*  LEAVING  THt  MAIN  BANK 

PH  12  =2  .0  *MSS  *  *0  .4  *C'ST/  KST 

TTT=<  T«S*(  PH  I  2- 1.01*2.  0*7YY)/  I  PH  12*1  .Oi 

*  C  C  M  FU T  E  TFE  TEMFEPATIJPE  OF  THE  FLUE  GAS  LEAVING- 

*  THE  ECONOMIZER 

PH13*2.C*CTL/nTL 

TUU=  ITTT  *1  °H1  3- 1.0*2.  C*TX/J  /{PHI  3*  1.0) 

*  C  3  M  PL  TE  ECCNCWI  IEF  ENERGY  TRANSFER 

*  FLUE  «AS  T)  TOPE  METAL 

Q  7*M  TT*C TL* ( TTT-TUU) 

») 


CS  M007  10 
CSMCC720 
CS'*  0o7  30 
CS  «CC740 
CSMCC750 
C  SM  00760 
C  S  MGG7  70 
CSM00780 
CS  MC0790 
CSMCOBOO 
CSM0C810 
CSM0U320 
CSMCC630 
CSM00840 
CSMCG850 

■:  sm  cc66o 

CS000870 
C  S*C0880 
C  SM 00890 
CS  M00900 
CSM0CS10 
CS  M  CC920 
CS  MO 09  30 
CSMCCS40 
CSMQU950 
CS  VC0960 
C  Sm  C  CS 70 
CSM0U960 
C  $  *00990 
CSMC1000 
CSMC1010 
CSMC1020 
CSMC1C30 
CS“01040 
C  S  y  C 1 C  50 
CSMJ  1060 
CS  *31070 

csMciceo 

CS'JC  1090 
C 5*011 00 
C  SM  c  1 1  1  0 
05  '•*  J  1  120 

CSm:ii30 

CSMC  1  140 
C j "0  1 1 50 
C5*C1UC 
CSMJ117G 

cs  Mcueo 

CSMCllSO 
CSM J1200 
C S  *  0 1 2  10 
CGV0U20 
C*  ”01230 
C  5*C1 2  40 
CS”  C  12  50 
CS  * J  12  60 
C  S  *  C  1 2  70 
CS"ciiao 
CG  *01  2 90 
3 S” c 1 3  CO 
CSM-.UJ10 
C  r  *01320 
C  Su  0  l  3  20 
CS  *01340 

i;  s  m  c  i  ’ ;  o 

OSMQ  12tO 
CS  MC 1 3  70 
C  S>'  C 1 3  EC 
CS  M  0  1 390 
CS*01400 
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FILE:  CS*° 


FORTRAN  PI 


NAVAL  POST 6ft  ACUAT  E  SChOOL 


* 

* 

* 


LSE  CONTINUITY  RELATIONSHIP  TO  COMPUTE  THE  MASS 
FLOW  -  ATE  UP  TH  E  5  ISE°  S 


VU.=MJJ-ORHCJL*VOLJL 
►  Ff=MCC-CPHl?CF*VnLCF 
ORAL  E  =  ufc  ft  1 VtJP'l O.fthOJU 
U"  42  E=  CE~  IV(  OR  A20.P  HOOF  ) 

FRCCECLRt  CRFCJL  =  MLT?1  IDRA14) 

IFIK.EEP.NE  .1  IOC  TO  1 

0PF3JL=  (  03  4  13*03  i  17034  160“  \  1  5+ OR  A1  40J  A  l  303  A  1203  4  l  1+  .  .  . 

GftA31*CiftAi2  03A33  03  434  OP  435  03  A36  ♦U<s  A3  7  +0=  A  3d  I  /  16. 0 
Q345j=DR452 
03452=  Cft  A51 
0R451=C3A31 
C3  A3  1=03  A  32 
03  A3 2= C3A33 
DR  A3  3=0R  434 
OR  A3 4=  CP  A35 
U  3  A3  5=  OR  A36 
OR  A3 1=03  A  3  7 
OR  A3  7=  OR  A38 
Oft  A3  £=ORAl 1 
CRA11=  CR  A 12 
0RA12=C3A13 
03  A  1  3=03  A  l  4 
03  A14=  C3  A  15 
JkAI  5  =  0.RA16 
ORA  lc=OR A  1 7 
URA1 7=  OR  A1 8 

1  CCNTINLE 
t  'JCPRnC  EDU3  F 

FPCCE  CURE  CPFCCF=FIL'"R2<JFA24) 

IFIKFEP.N'E.  II  SO  TO  2 

CRFDCF=(  C«A2oOR  427  0=  A2oO°  A  2 50*  A  24+03  A  23+0  ft A 220ft A2 1+ . . . 

0RA41  +  CRA42O3A4  3  *03  449  +0PA45  +0  3A4u+DPA4  7  +  03A4d )/ 16.0 
03  463=03*62 
CRA6  2=CRA6l 
0RA61=0RA41 
OR  A 4  1=  OR  A 42 
ORA42  =  0RA43 
03  A 4 3=03 A  44 
CP  A4 A=  LR  A45 
Oft  A4  5  =  0ft A46 
OR  A4 6s  C3  447 
CKA4  7=  03  A4  8 
OR  A4  E=DR A2 1 
OR  A2  1=  CP  422 
0RA22=0RA23 
03423= 0RA24 
CftA2  4=CRA25 
OR A2  5=0R  42  6 
C°  426=03427 
OPA27=CRA2a 

2  C  ON  T I  N  LE 
E'J  C°PQCtCURE 
* 

*  COMPUTE  THE  average  GENS  i  t  y  in  tfe  risers 

« 


RHOJL=(  1.0/L jl)  =  I (LKL/I ALL* VFG) )* ALOGU  X L L / VF I *VFG  +  1 .01... 
♦  FHCJJ*L  M ) 

RHCGF  =  (1 .O/LCFI* ( (LEFA 1XFF*VFG)1 *ALOG( (XFF/VF  I  *V FG ♦ 1  .0  I  .  .  . 

♦PH1CC*L JE I 
PHCFF*  =ECf-XFMPhOFG 
»MOLL=RH 0F-XlL*3hqF0 

•  '•tmpuTE  THE  MAI  1  RAN<  ENcRGY  TRANSFER  - 

•  TUBE  HETAL  TC  WAIN  BANK  M  IxTU5  E 

• 

C6  =  K  Y  *  c  S  AT**(4.0/ 3.C)*<  TYY-TS4T) **3. 0 


C  SPC141C 
CSM  0 14  20 
C  SRC1430 
C5MC1440 
CS "0  L450 
C  S  MC 14  60 
: 3" 014 70 
C3  v  J  1 4  80 
CSPC14S0 
CSOliCO 
CS  “01  510 
; SM  Cl  520 
CSM0I530 
CSMQI540 
CSO1550 
CSO1560 
CSMC1570 
CS“0  1530 
CS VC1590 
C S“ C 1 6C0 
C3 “0 16  10 
CS  “C 1620 
C  S-  C163C 
CS  R 0 1640 
C  S  *■'  C  1 6  50 
CS  C  It  60 
CS* 01670 
;  S'-  C1680 
f.S“  J1690 
C  3 “Cl  7  00 
C  SO  17  10 
CS  >'71720 
CSMC1730 
C3  '■*  0 1  740 
C  3 PCI  750 
:s"ci7to 
C  3  V  J 1  7  70 
CS-CI 780 
C3“0  1790 
CS  01800 
CS“C  1  6  10 
CS''-C1820 
:s 01830 
:  s 4 £ ia40 
CS  016  50 
C SVC  i860. 
C  J  '01870 
CS  01.180 

c  S’4 a  e so 

CS“  01 900 
C  S>'C19  10 
C SVC  15  70 
CS  O  19  39 
C  SOI SAO 
CS" J 1950 
CS VO  i960 
CS-C1970 
CS“^1980 
CS “J 1990 
C  S'4  C2CCC 
C  >*'02  0  10 
C  S yC2  020 
SMC2C30 
C3O2J40 
C  SVC.’OSO 
C  SO  20  60 
C  3  Oc  J  70 
csMicceo 

CS  ’’  C2090 
-SOilOO 
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FILE:  CSMP  FORTRAN  ®l  NAVAL  PC3TGR/CLATE  SCHOOL 


*  COMPUTE  THE  SCREEN  RISER  ENERGY  TRANSFER  - 

*  T'JEE  METAL  TO  SCREEN  RISER  MIXTURE 

* 

02=KV*PSAT**( A.C/3.C)*! 7  V  V-TSATI ** 3. C 

• 

*  COMPUTE  THE  RISER  AVERAGE  ENTHALPYS 

x 

F0F=H0CH-XFF*HFG/2.0 

1-  JL=HJ  J+XLL*HFG/2.  0 

* 

*  CCMFLTE  THE  RISER  OUTLET  ENTHALPYS 

* 

hf  F=FF+X  FF*FFC- 
F  L l=HF ♦ XLL*HFG 

*  C3MFUTE  THE  RISER  EFFECTIVE  HEIGHTS 

* 

iC  E=  Z  CF-L  E  F 

2- i>  =  Z2l-lKl. 

ZNL=ZJL-ZJK 
ZEF=  ZCF-ZCE 

♦ 

*  COMPUTE  THE  SPECIFIC  VOLUME  OF  THE  ECONOMIZER 

*  C'JT  IE"  LIQUIO 

* 

V  3  8=  .0160C46o-.C0CJ020  146*T38-*-.  0Q0C0C03651  l*TPR**2.0... 

-R.142E-11  *TE3**3  .)+i..40i£-l.J*T3H**4  .0-1 . 143E-  16 *T 33 **5. 0 .  .  . 

*  c.  0  24E-  2C*T33**o.O 

X 

*  CCMRLTE  the  RAtE  tO HAT  I CN  FOR  RISER  OUTLET  QUALITY 

* 

ChX  LL=  (MCGMFJJ-fF-xLL  xfFG/2.0  )+Q tr- MLL *XLL  *HF  G /2. 0>  /  l RHO JL* VOLKL ) 

HXLL  =  I NToRL (HXLLC  .CFXLL  I 

aLL=<H>LL-FF)*<.0/3FG 

CF/FF=(MCC*<  FCC-FF-XFF*FFG/2.0  )+  02-M  FF*X  FF*  HFG  /  2 .0  )  /  I  R  HOOF  *  VOLFF  ) 
UXFF  =!  MGRLHXrFO.OHXFF) 

XFF= (FXFF-FF  1*2  .J/FFG 


*  CCMPLTt  ThE  RISER  BOILING  VOLUME 

« 

VCLEF  =  VCLC  F*LEF/LCF 
VOLKL=VJL  wLX'lKL/LJL 

4 

*  CCYF'ITE  THt  NON'EOIL  IMG  LENGTH  OF  THE  RISERS 

* 

10E  =  LCF*  1 1- F-H  C C  1  /  <  (HF+aFF*FFGJ-HDO  ) 

U*=LJL*  IHf-rUJl  /  ((fF-*aLl*EFG)-HJJ  I 


COMPUTE  The  RULING  LENGTH  op  the  risers 


* 

* 

L<  L=  L JL-L JK 
l£F=  LCF-LCE 

* 

X 

* 

* 

*  COMPUTE  THE  MASS  RATE  EQUATION  FOR  STEAM  CONDENSING  IN 

*  ThE  CRUM 

* 

MCCNC=560.'73*(°''M/{TMM*4o0.0)**0.5-0SAT/(TSiTF<i60.J)**0.5)... 

♦  . C2568 

* 

* 

* 

COMAS  L=MLL  *  (  l.J-XLL  )  +MFF*(  1  .C-  XFF  J  FMCOND  M  B3 0- MCC-  MSG 


CCMF’IT  £  THE  RATE  EOUATIHJ  FOR  DRUM  LIQJID  MISS 


COMPUTE  THE  DRUM  LIQUID  MASS 


CSMC21 10 
CSMQ2120 
C  S^ 02130 
c S' 02140 
CS  v02 1 50 
C  SM021 tO 
CSMC2170 
CS ”02180 
C  SM  C  2 1  SO 
CSM02200 
C  3M  022  10 
C  SMC2220 
CSM02230 
C  S  *02240 
C  SM02250 
CS  *02260 
CSMC2270 
CS  M  C22  80 
CSMC2290 
C SMC23C0 
CSM32310 
C  SVC2  320 
CSMC2330 
CS "02340 
C  S  MO  23  50 
C  S  m  C  2  3  60 
CS  *023  70 
CSMC2380 
CS  " C23  SO 
CS  *02400 
:  s'‘  C2sic 
CS  *  02420 
C  S *  C2430 
CSMC2440 
CS  *02450 
CS*C246C 
CS"  0 24  70 
C  S  M  C24  30 
C  S" L24S0 
CS  *'02500 
C  S  y 025  10 
C  S’M  0  2 520 
CS  '*02530 
C  $  v  C  2  5  4  0 
C5MC2550 
CS* 02 560 
C  SMC25/C 
CS'*  02  5  80 
C  S*C25S0 

:  s*-*  c2tco 
CS  MO 26  10 
: S  *02620 
!_■  0  26  20 
CS  M 02 040 
CS*C2650 
CS " C26  60 
CS '026 70 
C  S  *  C  2  6  6  C 
CS  M  0  26S0 
CS  *02  700 
SS'*C2710 
CS' 02720 
CS'*  02730 
C  S  *•*  0  2  7  40 
CS  w 02  750 
C  S  *  C  2  7 1 0 
CS" 02  7  70 
C  S  MO  27  30 
C  S'*  C2  7S0 
csuu28oa 


8.? 


FCRTRAN  PI 


NAVA l  posTGRADlaTE  SCHOOL 


FILfci  CSHP 
j«lASL  =  lNTGRL<DVA$LO  .DDMASL) 

J  CD“P'JTE  thH  S  AT  6  E  00  ATI  ON  FO5  EN  E3  GY  IN  T  hE  J3U“  LIQUID 

CJ^OnL  ="  LL  *  (  ].J-XLL)*HF*'IFF«(  l.O-XFF  )  =*HF  *MCOND*HFG*M34*H8;3... 

—  i»CC*hCC'~M'>»*Hvj>j 

*  COMPUTE  THE  ENERGY  IN  THE  DRUM  LIQUID 
C“CHL=  INTG»L  (  DMDHLO.DDMDHL  I 

*  COMPUTE  THE  ENTHALPY  CF  THE  DRUM  LIQUID 

* 

OH=DMDHL/CVASL 

* 

# 

*  COMPUTE  THE  DRUM  SPECIFIC  VOLUME 

• 

CRYUNO  =  PCU*  (Mll  +  VFF) 

PROCEDURE  RI  3c2=F  HTR4  (CRYUND) 
p  1S£1=GEL  1Y{  <aU,K.  ISTIM.CR  YUND) 

R I SE2*CRYU MD 

I F  (  VALVE. GT..  51IGQ  TO  55 
R  IS  e  =  R  IS  E2 
i>C  TC  57 

55  P  1 S  t  =  R  IS  E  1 

57  C  CM  11  NUE 

ENC°’nC  FCiJR  t 

CC  MOV  =  (('/FF-*VLL-’S4J*VF*,CI*(  mFF«-mll  )  *V  V  *MBg  Q«  v*3- 1  MCC+  M  GG)  *  VF  «-MC  OSD 
*  VpG-Pl  Sc*  VV) 

CMCV*INTGRL<0  IDVO.DOMDV) 


COMPUTE  THE  OR'JM  LEVEL 

.EVEL*(D»CV-VCL,.RY/2.J  j/  < LST “OM*CSTMCm ) 

CO^UTF  THc  COWNCOMER  ENTHALPY 


HCO=  1  (  FFF*PLL-M83)*HF»ME3*F33  if  l  MFFmMLL) 

FGF=  FC  D 

FCC=FCC 


FCC=HCO 


FOG3  HGH 


* 


*  COMPUTE  THE  COwMCCMER  SPECIFIC  VDL’J“E  AND  DENSITY 


♦ 

V"D  =  (  (  PFFMMLL-MBa)  «VFt  MeS*  V3  0)  /  (  MFFmmll.' 

=  hicc=  wo/vro 

RFCGH=RHCCO 

* 


*  COMPUTE  THE  SATURATION  PRESSU3E  A  NO  TEMPERATURE 

*  CCRFtS  FGNCINU  TC  THE  CTWNC'mER  ETHALPY 


PS  AT=£xP  ( (  /LOG  (FC  at  )-<*  .4>7C3  )/ .  204  52)  . 
TSAT=E  >P  l  1.22151*  AlG»<  PSAT  J  +4.77123  H 


«  CCMFUTf  The  £  |T  FALPY  CF  TFE  LIQUID  ffJTERING  TFE 

*  MAIN  BANK  R1  SER 


H JJ=HGH+C9/MH 

* 


*  COMPUTE  ThE  RATE  EQIIATICN  FOR  TFE  MAIN  SANK  AND 

*  SCREEN  RISE*.  rUeE  MET AL  TEMPERATURES 

* 

DT VV  * { Cl -C2J/ <  FASSV*CFVI 
DTYY=(Q5-QtJ  /(  MASJY’CPYJ 

» 

*  CCMPLTE  TFE  SCREEN  AND  “AIN  RANK  RISER  TUBE 

*  ME  7 Au  TEMPER  A  rURE  S 

* 


CSMC281D 
C  Sv  02820 
CS “02330 
CSyC2840 
CSMC2E50 
CS  ‘••C23  60 
CS“C28  70 
CS  M  02  3  SO 
CSMC2S90 
S  M  C  2  S  C  0 
C  >“029  10 
C  SMC2S20 
C3“C29  20 
CSMC2940 
C  S*’  C  2S  50 
CSV 02960 
CSMC29  70 
C  S- 02980 
CS “02990 
C  S “C  2C  CO 
CSVC 30  10 
CSMC3020 
C  S“ C  a  C30 
CS“ 03040 
C  S “03 050 
CSMC3C6C 
CS “C30  70 

:  sm  C3  c  80 

CS “03090 
; S“C3ioo 
C  3  ‘ '  C  2  l  10 
CSMJ3120 
C  S “ C  J 1 30 
CS-03140 
C 3 “Ca l 50 
C SMC  a  160 
C  S v 03 1 70 
C  S “C  3 1 SO 
:  S“ 0  3  190 
CS “03200 
C  S  C  2  2  1 0 
CSMC2220 
CS  mo 32 30 
C  SM  C  32  4C 
CS " 0  2 2  50 
C  SF  Op260 
C  >-0  32  70 
CS  VJ3280 
C  SM  C  3  2  SO 
CSU J  3 300 
CS “ 0>3  1 0 


C  S  '  c J  3  20 
CS- J j 3  30 
C  S  “  0  j  3  40 
:  S’1  0  3  3  50 
CS  yOi360 
C  SM  C  2  2  70 
CS’-’  03380 
CSM; JA90 
C  S“ C34C0 
C >“03410 
C  S-C2420 
CSM  C  34  20 
CS“C  3440 
C S“ 1 34  50 
r.s“o>460 
U  S“C34  70 
CS“o;4P0 
CS “0  1490 
C  S “C35CC 


FORTRAN  R  1 


NAVAL  POSTGRADUATE  SCHCOl 


r 


ULf-  CS-P 

rvv=INTGRL<  TVVC.OTVV) 
-  YY*  1 NTG  BL  (TYYO, CTYY) 


COMPUTE  THE  HEAT  TRANSFER  FOR  THE  DESUPERHEATER 
-STEA*  TO  TUBE  NETAL 


09=CNP  TNN-TPP ) 


« 

* 

* 

TPP=  (  TNN-T ZZ  i  /  ( E  A P (  KNP  /( CN°  *MNN**0.2  ) )  I  TZZ 

* 

* 

* 

* 

TGOTSST 
1HH=TSAT 
TC  C=  TS  AT 
TOC=TS  AT 

* 

* 

* 

* 

09 9=  K  Z  *  ( TZZ-THH) 

* 

* 

* 

* 

CT  ZZ=  (09-0  99  I  /{  MA  SSZ  *C  PZ ) 

* 

* 

* 

TZZ  =  INTGRUTZZO,CTZZ) 

* 

*  compute  THE  ECONOMIZER  ENERGY  TRANSFER-TUBE  METAL 

*  TO  FEEL  WATER 

* 

CS=MAA*CAE*(TEc-TAAI • 


CCVFLTE  THE  TEMPERATURE  OF  THE  STEAK  LEAVING  ThE 
CtSUP  EP  HEAT  ER 


SET  THE  DOWNCOMER  ENTRANCE  AND  EXIT  TEMPERATURES 
EQUAL  TO  ThE  SATURATION  TEMPERATURE  CORRESPONDING 
TC  CRUK  ENTHALPY 


COMPUTE  THE  ENERGY  TRANSFER  FOP  THE  DESUPERHEATER 
T'JEE  METAL  TC  DRUM  LIQUID 


COMPUTE  THE  RATE  EOUATI  ON  FOR  DESUPERHEATER  TUEE  METAL 
TEMPERATURE 


COMPUTE  THE  DES  LPEKHE AT  ER  TUBE  METAL  TEMPERATURE 


CDMPLTF  THE  FEED  TEMPERATURE  AT  CUTLET 
CF  ECONOMIZER 


» 

* 

* 

T8B*  (TAA-TXX)  /  <SXP(  <X/  (CA°  «M.\;  **c.  2)  ) )  TXX 
* 

*  COMPUI'E  THE  RATE  EQUAT1CN  FOR  THE  ECONOMIZER  TLEE 

»  R£T  A  L  TtMPEPArUFE 

* 

OTXX  =  (  08-0  71/1  MASSX*CPX) 

4c 

*  COMPUTE  THE  ECONOMIZER  TUEE  METAL  TEMPERATURE 

* 

T>>=I NT^PL (TXXC , CTXX  J 

4 

«  CCMFU-E  THE  SUPERHEAT  EC  STEAM  OUTLET  TEMPERATLRE 

• 

TNN=  (Tmm_T»(W  J  /  (  EX3(  <W/  (CMN  4MNN  **C.  21  JJ  4-TWW 

* 

*  COMPUTE  THE  SUPERHEATER  ENERGY  TRANSFER-TUBE  "£T  AL 

*  TO  ST  cAM 
C4*C  KN*KKK* ( INN-TMM) 

* 

* 

* 

* 


C0M  pip-  g  THE  RATE  EQUATION  FOR  SUPEPHEATER  TUBE  METAL 
T=  MPE  RAtUR£ 


CT  W W  *  ( Ci -04  ) /  (  KAS3W*CFW) 


4 

* 


COMPUTE  THE  SU°  ERH  EA  TER  TUBE  METAL  TEMPERATURE 


CSMC3510 
CSM03520 
CSMCJ5  20 
CS  M03540 
CSMC3550 
CSMC3560 
CSM035  70 
CSMC35E0 
CS’*  C  3590 
C  SMC36 00 
CSMC3610 
C  S  MQ3620 
CSM0363C 
C  Su  0  36  40 
CS«03650 
CSMC3660 
CS MO 36 70 
CS  "03680 
C  SM  C36S0 
CSM03700 
CSM03710 
C  SM  0  3  720 
CS "03730 
C  SMC3740 
CSm  03750 
CS  "03  760 
C  SM  C  3  7  70 
CSMC37EO 
CSM03790 
C  SMC3EC0 
CS "03810 
CS“C2820 
CS'1  02  6  30 
CS  u03840 
CSMC3E50 
CS  M  03  860 
CSMC3870 
CSM  C2BE0 
CS “03890 
U  S  "0  39  00 
C  S  M  039  10 
CS "03920 
CSMCJ930 
CS " C3940 
CS " 03950 
CSM  02960 
CS  "0  3 9  70 
C  SM  C  39  BO 

CSMC2SS0 

CS" 04000 
C  SM  C  4  C  1 0 
CS  “  J  40  20 
CSM04030 
CSMC4040 
CS  04050 
C  S"C4CfcO 
CS'-'  040  70 
CS  "C-»080 
C  SM  C4090 
CS  04  1  00 
C  S "C4 1 10 
CSM0il20 
CS "04130 
CS"C4140 
CSMC4150 
CSM041&0 
CSUC417C 
CS"04l  tiO 
CS "C ♦ l 90 
CSMC4200 
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ULE:  CS“P  FORTRAN  »1 


T*.*=INTGRU'UVC,OTWI.) 


CCMFUTE  Th£  IMPLICIT 
FLC*  RATE  FRCM  E31LER 


EOUATIO'J  FOP  TOTAL  STEAM 


l-ff-1  £FL1*“M0  "MlfPI 

*  Cn  “  FUT £  THE  SJPEPHEATEP  OUTLET  PRESSU°E 

3*'N=SO°T(DSAT  ♦  *  2  •  0~  2.  0*KONi*KaNl*MMM  4*2. 0) 

““MI  II  *PAMKCI\4»VALVE 
mm«j  j  =  M\fJO 


*  EOlATE  THE  FLUE  GAS  MASS  FLOW  RATES 

* 

“PRaMQQ 

••■SS^MPP 

►  T  T* “S  S 
MLL*MTT 


SET  kIScR  INLET  CEOS  ITY  EQUAL  TO 
SATURATED  LIQUID  DENSITY 


PH  CEE  =  RH  CF 
5hOJ  j=RHCF 
PHCGC.PhOEE 
RHCKI>  =  fiHC  JJ 


* 

* 

»HH*MGG 

‘CC=“0C 

“G  \j  —  “  J  J 

MOOaMFFC 
“  J  J*  “LLO 

« 

* 

* 


ECUATg  CCWNCCMtP  FLrW  RATES  T"1  °ISE=  ENTRANCE  FLOW 
RATES  A r>.J  SET  RISER  ENTRANCE  FLCW  S AT £S  TJ  THE 
FLOw  RATES  AT  THE  INITIAL  TIME 


COMPUTE  THE  DERIVATIVE  OF  AVERAGE  RISER  DENSITY 


4 

* 


COMPUTE  THE  STEAM  MASS  RATE  EQUATION  FOR  TMt  STEAM 
DR1  IF 


ODSTms  *FF4mFF+-XLl*“LL-  “CCNO-“MM 


* 


*  COMPUTE  Thfc  CRIJM  STEAM  MASS 

4 


CSTM*INT  CP l  t  CS  T  “0  *  DOST  M J 

4 


*  COMPUTE  THE  VOLUME  OF  STEAM  IN  THE  STEAM  ORJM 

* 


VnLST“=VQLCRM-OMASL/RHOF 

* 

*  COMPUTE  THE  DENSITY  OF  STEAM  IN  THE  STEAM  C3'JM 

4 

RHCST“=CcT“/VCLSTM 
PHOmMsRHCSTM 
VV“M  =  1  .0/PHCMM 
VFGMK*VVF8-VF 

* 

*  COMPUTE  THE  STEAM  0=  UM  STEAM  OUTLET  PRESSURE 

* 


noESS“=524.C/(VFGMM*.  1) 

PRCCEOURE  FMM  =  F1LTP5  IPRESSM) 

I  r  |  o*F  SSm.LT.  J.  OJCALL  UE3US  (3*0.0) 

FMVx  PRESS'* 

t  NJFPCCfcCLRE 

•  COMPUTE  THE  STEAM  OR  JM  STEA'-*  OUTLET  TEMPERATURE 


C5M  C42  10 
CS  “Cm220 
C  S“  C42  30 
CS  m  c**240 
CS“G42  50 
c  suC‘.2eo 
CS “042  70 
CSMC4280 
CS  u  C42  SO 
CS“0**300 
C  S’M  C4310 
CSM 04320 
CSMC4330 
CSM04340 
CS “04350 
C  SMC4360 
CSMQ4370 
CS  “04380 
CSMC42S0 
CSMC44C0 
CS  “04410 
C  SM  C442C 
CS‘*  04430 
CS“C4440 
CSMC4450 
CS“04460 
CS “044  70 
CSMC4480 
CS  “04490 
:s,jc45Co 
CSMC4510 
CS“C4520 
C  SM  C4530 
CSM04540 
C  S“C4550 
C  S‘’C4560 
CSM04570 
C$“04580 
C  SMC45S0 
CS“04600 
CS''C4tlO 
CS "04620 
CS “C4630 
C  S-"C4fc4C 
CS  v  04650 
CS“C4b60 
CS“C4t70 
CS  “04680 
CS“C46S0 
CS“ 0*700 
C  S“G-*7 1 0 
CSMC4720 
C  S “047  30 
C  S“0-.74O 
C  S’1  04  7 50 
CS “04760 
C  S“ 04 7  70 
CS C4  7 EO 


CS “04790 
r  S“ C4  8C0 
CS  M  Cm  810 
CSMC4820 
C  SM04830 
CS “04840 
C  S“ C4  E 50 
CS  M 04860 
C  S“C48  70 
C  S“ 04 8  80 
CS “04690 
CSMC4900 
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FILE:  CSMP 


FORTRAN  PI 


NAVAL  POSTGRADUATE  SCHCOl 


TMM=£X°( . 22 15 l*ALOG( P'AMJ  *4.77123  J 

SOL  \c  FOR  THE  ECCNO*I 2ER  FEED  CUTLET  P'lT H  ALpY 
BE*IF5B*HAAO*C8)/M88 

FEECRATE  EQUATION 

8B*M880 

Two  Phage  FLOW  "UlTI°L1ER  E3U AT  IONS 

F6FAVE=24.704*XFF**2  »J  —  6  • 5  0  Lb  *X  F  F*  .  9  77t> 

R  jS A VK=24. 754*  XL L**  2.0- t. 50  i  6*  X L  L*  .  9776 
;aat=l3.‘i5tt**Ff*»2.»J>lo.  49 44* A FF~. CO 007 
RACL  J  =  15. 4 56-*  *XLL**2  .  J  *13  .'4 44*  *A  LL - .30  30 7 
"FRl  CE=-  34.  Ct22*xFF**2. 0  2  3.  71  6  4*XFF  +.8  734 
FFiCK  =  -J‘t.Oo22*XLL  .C»2 3 . 7 1 6 4*XLL ♦. b 7 34 

STATE  POINT  EOLATIONS 

SAT=*DH 
r=HSAT 

hFL=922.15-0  .=»35l6*PSAT4i..7i7£-04*oS  AT**2  . C-4 . 2 19E-0 £.  .  . 

*PSAT**3. 0 
HV  =  HSiT*hFG 

FriCF=63.8-0.31781*TSAT  +1  .1 32  E-05  *TS  AT**2  .0-6 . 786E-03  .  .  . 


CSMC4S10 
CSM04920 
CS*C4930 
CS"  C4940 
CS  "04 950 
C  SM  C49t0- 
CS"C4970 
CS  "049  80 • 
CSMC4S50 
CS-’OSOOO 
CS*C5010 
C  S"  C  5C 20 
CS "05030 
C  S  " C6C40 
CS  ‘"050  50 
CS  *05060 
C  SMC 5 070 
CS "05080 
CS*05090 
C  SMC 5  100 
CS"C5110 
C  SM35 120 
CSMC5I30 
CS"D5140 
CSMC5I50- 
CS ’<05160 
CS  *05170 
CSMC51E0 
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G;  Anticipated  Performance 


APPENDIX  C 


Design 


Endurance 

Rated 

Full 

Maximum 

Intermittent 

Endurance  Boiler 

2  Boilers 

Power 

Power 

1  Boiler 

Overload 

Rate  of  Operation  -  Per  Cent 

56 

83 

97 

114 

120 

Total  Steam  Generated  ib/hr 

195,000 

290.000 

340,000 

400.000 

420.000 

Superheated  Steam  Ibs/hr. 

155,000 

235,000 

265.000 

300,000 

320.000 

Desuperheated  Steam  ibs/hr. 

40,000 

55,000 

75,000 

100,000 

100.000 

Boiler  Drum  Pressure  psig 

690 

690 

690 

690 

690 

Superheater  Outlet  Pressure  psig 

665 

630 

610 

580 

570 

Superheater  Outlet  Temperature  °F 

912 

914 

980 

898 

895 

Desuperheated  Steam  Outlet  Pressure  psig 

651 

606 

558 

480 

470 

Desuperheated  Steam  Outlet  Temperature  °F 

635 

659 

680 

690 

690 

Economizer  Iniet  Pressure  psig 

715 

726 

733 

749 

754 

Economizer  Inlet  Temperature  °F 

2S0 

250 

250 

250 

250 

Economizer  Outlet  Pressure  psig 

703 

705 

703 

708 

710 

Economizer  Outlet  Temperature  °F 

355 

372 

372 

377 

380 

Casing  Air  Inlet  Temperature  °F 

100 

100 

100 

100 

100 

Total  Air  Flow  lb/hr 

242,229 

368.092 

433,589 

513.154 

540.952 

Total  Oil  Flow  Ib/hr 

14333 

21 ,781 

25,656 

30,364 

32,009 

Anticipated  Efficiency  7a 

85.7 

84.4 

83.8 

83.0 

82.7 

Guaranteed  Efficiency  % 

85.7 

84.4 

83.8 

83.0 

82.7 

Rad;ition  and  Unaccounted  for  Losses  % 

.99 

.99 

.95 

1.01 

.98 

Excess  Air  % 

15.0 

15.0 

15.0 

15.0 

15.0 

Carbon  Dioxide  ~c 

13.0 

13.0 

13.0 

13.0 

13.0 

Number  of  Burners  in  Operation 

6 

6 

6 

6 

6 

Throttle  or  Non-Throttle  of  Air  Doors 

N.T. 

N.T. 

N.T. 

N.T. 

N.T. 

Draft  Loss  -  Total  Inches  Water 

13.74 

33.98* 

47.96 

69.51 

76.22**  93.48 

Through  Double  Casing 

133 

3.09 

4.28 

6.00 

6.67  10.22 

Through  Burner  Register 

2.60 

6.00 

8.50 

11.50 

13.00  15. SS 

Through  Boiler  &  Superheater 

4.83 

12.29 

18.14 

25.61 

28.15  38.32 

Through  Economizer 

4.98 

12.60 

17.04 

26.40 

28.40  29.06 

Gas  Temperature  Leaving  Superheater  Screen  °F 

2471 

2594 

2642 

2686 

2699 

Gas  Temperature  Leaving  Superheater  °F 

1826 

1961 

2021 

2085 

2103 

Gas  Temperature  Leaving  Main  Bank  °F  • 

694 

774 

812 

856 

872 

Gas  Temperature  Leaving  Economizer  °F 

373 

427 

453 

484 

498 

Heat  Release  KB/H./Sq.  Ft.  Radiant  Heat  Absorbing 

Surface 

447 

680 

800.6 

948 

999 

Heat  Release  KB/Hr/Sq.  Ft.  Total  Heating  Surface 

15.5 

23.6 

28.0 

32.9 

34.7 

Heat  Release  KB/Hr/Cu.  Ft.  Furnace  Volume 

183.9 

279.5 

329.2 

389.6 

410.8 

Furnace  Heal  Absorption  KB/IIr./Sc;.  Ft 

128.4 

171.0 

191.1 

213.19 

221.5 

Heat  Absorption  First  Water  Scicen  Row  KB/Hr/Sq. 

210 

231 

255.3 

288 

256 

Ft.  (Max.) 

Heal  Absorption  Maximum  KB/I!r./Sq.  Ft 

210 

231 

255.3 

288 

256 

(Furnace  Screen) 

•Draft  Losses  (Full  Power)  Based  On  241  Cu.  Ft.  of  100°F  Air/Lb.  of  Fuel  Oil 

••Draft  Losses  (Overload)  Based  Oil  241  Cu.  Ft.  of  IOC°h  Air/Lb.  of  l  ucl  Oil 

•••Draft  Losses  (Overload)  Based  on  2o0  Cu.  Ft.  of  68°F  Air/Lb.  of  Fuel  Oil 
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H.  Anticipated  Metal  Temperature  Degree  F. 


Endurance 

2  Boilers 

Water  Screen  Tubes  Outside 

644 

Water  Screen  Tubes  Inside 

554 

Superheater  Tubes  Outside  -  Maximum 

1022 

Superheater  Tubes  Inside  -  Maximum 

1044 

Maximum  Inner  Casing  Temperature 

SOO 

Maximum  Outer  Casing  Temperature 

at  way  of  Structural  ties 

350 

Maximum  Outer  Casing  Temperature 

145 

Design 

Rated 

Maximum 

Intermittent 

Endurance 

Boilei 

Full  Power 

Power 

1  Boiler 

Overload 

656 

661 

667 

680 

554 

554 

554 

555 

1046 

1049 

1048 

1048 

1020 

1019 

1014 

1012 

825 

835 

850 

870 

350 

350 

350 

350 

145 

145 

145 

145 

Tube  Data 

OX>. 

M.W.T. 

No. 

Material 

Side  Wall  and  Roof 

2” 

.134” 

71 

MIL-T-16286  CL.  A 

Rear  Wall 

2” 

.134” 

54 

MIL-T-16286  CL.  A 

Front  Wall 

2” 

.134” 

22 

MIL-T-16286  CL.  A 

Screen  Bank 

2” 

.134” 

102 

MIL-T-16286  CL.  A 

Superheater 

1.5” 

.120 

22 

MIL-T-16286  CL.  L 

Main  Bank 

2” 

.134” 

34 

MIL-T-16286  CL.  .* 

Main  Bank 

1” 

!/■> 

00 

p 

2808 

MIL-T-16286  CL.  A 

Economizer 

2” 

.180” 

182 

MIL-T-16286  Cl  .  a 

Desuperheater  (in  Water  Drum) 

2” 

.220” 

6 

SA-268  TP430 

Downcomers 

• 

8  5/8” 

.483” 

1 

Schedule  80  Pipe 
ASME  SA-106-B 

Downcomers 

10  3/4" 

.519” 

5 

Schedule  80  Pipe 
ASME  SA- 106-15 

Downcomers 

12  3/4” 

.601 

Schedule  80  Pipe 
ASME  SA-10n-B 

Risers 

6” 

.500” 

7 

ASME  SA-IOlvP. 

♦These  tubes  may  be  MIL-T-16286, Class  A  (Seamless)  01  MIL-T-1 7188  (Seamed)  electric  resistance  welded. 
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